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I.  INTRODUCTION 


A.  Statement  of  the  Problem 

Nearly  all  materials,  including  soil,  lose  strength  with  repeated 
loadings.  The  nuaber  of  repeated  loadings  required  to  reach  failure 
increases  with  decreasing  cyclic  stress  intensity.  When  a  dry.  loose 
sand  is  subjected  to  stresses  sufficient  to  cause  intergranular  slip  or 
particle  fracture,  the  rearrangeaent  of  grains  leads  to  volumetric 
coapaction.  Voluaetric  compaction  can  be  induced  by  a  single  large 
loading  or  a  series  of  loadings  of  smaller  intensity.  When  a 
saturated,  loose  sand  is  subjected  to  stresses  of  sufficient  intensity 
to  cause  intergranular  slip  or  grain  fraoture,  compaction  is  inhibited 
because  the  water  cannot  drain  instantaneously  to  accommodate  the 
potential  volume  change.  As  a  result,  the  porewater  pressure  increases 
since  the  water  must  support  the  gravity  loading.  The  larger  the 
increase  in  porewater  pressure,  the  greater  the  loss  in  strength  of  a 
oohesionless  soil.  Depending  on  drainage  conditions,  the  residual 
excess  porewater  pressure  can  persist  for  minutes,  hours  or  even  days 
until  enough  drainage  occurs  to  accommodate  the  volume  change. 

Porewater  pressure  increases  in  loose,  saturated  sands  can  be 
produced  by  sustained  low  frequency,  low  amplitude,  random  ground 
motions.  A  complete  loss  of  shear  strength,  as  a  consequence  of 
reduced  effective  stress  due  to  increased  porewater  pressures,  is 
termed  liquefaction.  This  behavior  can  lead  to  catastrophic 
consequences  including  landslides,  foundation  failures,  ground 


subsidence  end  eabankaent  fell ores.  Liquefaction  ss  s  result  of 
earthquake-induced  ground  notions  is  well  understood  end  bus  been 
doeunented  by  asny  researchers  in  both  the  field  and  laboratory. 
Liquef notion  froa  explosive-induced  ground  notions  is  not  veil 
docuaented  beyond  its  existence  and  the  aechanisn  is  not  well 
understood.  Field  dsta  are  often  United  in  extent  and  incoaplete  in 
detail  while  dynaaic  experiaental  laboratory  facilities  are  still  in  a 
developaental  stage. 

Kecognition  of  the  potential  effects  of  explosive-induced 
liquefaction  warrants  a  study  of  explosive-induced  liquefaction 
phenoaenon  and  related  effects  in  the  field.  Daaage  disproportionate 
to  the  aaount  of  explosive  used  or  of  a  nature  inconsistent  with 
previous  experience  nay  be  attributable  to  liquefaction.  For  exaaple, 
recent  re-exaainations  of  the  events  at  the  Paoific  Proving  Grounds, 
where  nuclear  explosives  were  detonated  in  the  1950' s,  seen  to  suggest 
that  liquefaction  any  be  the  priaary  factor  causing  the  unusually 
broad,  flat  crater  shapes  and  other  phenoaenon  observed  there.  Other 
ailitary  high  explosive  tests,  construction  blasting  and  snail  scale 
field  tests  have  also  shown  evidence  of  daaage  beyond  the  direct 
effeots  of  the  explosion. 

Field  observations  and  liaited  eapirical  relationships  are  the 
current  basis  for  evaluating  and  predicting  explosive-induced 
liquefaction.  However,  these  eapirical  relationships  do  not  have  the 
benefit  of  a  large  data  base  or  extensive  theoretical  foundation.  Of 
priae  iaportance  is  the  developaent  of  a  controlled,  systenatic 
experiaental  evaluation  of  liquefaction  in  both  the  laboratory  and 
field.  Such  a  study  would  vastly  iaprove  the  understanding  of  the 


The  primary  objective  of  this  investigation  was  to  experimentally 
and  systematically  evaluate  the  behavior  of  a  saturated  granular  soil 
subjected  to  dynamic,  one-dimensional,  confined,  compressive  loadings 
in  the  laboratory.  The  experiments  were  conducted  to  simulate  the 
field  loading  of  a  soil  element  located  near  the  detonation  point  of  an 
explosive.  Intense  compressive  loadings  having  millisecond  rise  times 
occur  in  this  region.  The  soil's  porewater  pressure  response  both 
during  and  after  the  passage  of  the  stress  wave  was  used  to  evaluate 
the  liquefaction  potential  of  the  soil. 

A  large  number  of  parameters  have  been  observed  to  affect  the  onset 
of  liquefaction.  Some  are  associated  with  the  soil  and  underlying 
bedrock  while  others  are  related  to  the  explosive  itself.  In  this 
study,  the  effect  of  variations  in  the  initial  relative  density  and  the 
initial  effective  stress  along  with  the  intensity  and  number  of  applied 
loadings  were  investigated.  A  projectile  was  fired  at  a  confined, 
water  saturated  sample  with  undrained  loading  conditions  to  simulate  an 
explosive  loading.  Experimental  measurements  included  the  applied 
loading  stress,  the  porewster  pressure  response  and  the  projectile 
impact  velocity. 

All  experimental  runs  were  conducted  on  saturated  samples  of 
Monterey  No.  0/30  sand  at  relative  densities  and  effective  stresses 
ranging  from  0  to  80  percent  and  86  to  690  KPa,  respectively.  Monterey 
sand  was  chosen  beoause  the  performance  of  this  material  under  static 


and  vary  low  frequency  (1  cycle  per  second)  loadings  is  well  documented 
in  the  literature. 

Analysis  of  the  data  included  an  evaluation  of  the  influence  of 
several  parameters  on  the  peak  and  long-term  porewater  pressure 
response  in  the  soil,  the  stress  wave  velocity  and  the  peak  particle 
velocity.  These  results  were  used  to  define  liquefaction  threshold 
limits  and  develop  empirical  relationships  for  predicting  porewater 
pressure  increases  in  saturated,  cohesionless  soils.  This 
investigation  has  lead  to  the  development  of  a  new  apparatus  for 
studying  the  dynamic  response  of  saturated  soils  and  has  provided  a  new 
approach  for  examining  compressions!  liquefaction.  This  study  has  also 
derived  and  documented  some  important  relationships  between  soil 
properties  and  compressions!  stress  wave  loading. 


II.  LITERATURE  REVIEW 


A.  The  Liquefaction  Phenomenon 

i.  Bffwtt  al  ItiflitliRjUga 

Observations  of  the  effects  of  liquefaction  have  been  made  in 
connection  with  earthquake  loadinga  of  soils.  The  1964  Alaskan 
earthquake  destroyed  harbor  facilities,  roadways,  buildings  and  homes 
in  Valdez.  Ala aka,  due  in  part  to  a  large  flow  of  liquefied  soil  (Seed. 
1968).  In  1964.  an  earthquake  in  Niigata.  Japan,  caused  extensive 
structural  daauge  aa  a  result  of  bearing  capacity  failures.  A  number 
of  buildings  sank  into  the  ground  while  a  sewage  treatment  tank  floated 
up  above  the  ground  and  an  apartment  building  rotated  intact  almost  80 
degrees  (Seed  and  Idriss,  1967).  Informative  reviews  of  earthquake- 
induced  liquefaction  phenomena  and  case  histories  are  given  by  Gilbert 
(1976),  Green  and  Ferguaon  (1971)  and  Seed  (1968,1979). 

Although  such  dramatic  and  widespread  damage  has  not  been  reported 
for  explosions,  it  is  not  unreasonable  to  expect  that  such  catastrophic 
consequenoea  as  those  associated  with  seismic  ground  motions  from 
earthquake  loadings  are  possible.  The  primary  sources  of  information 
in  this  area  are  from  indirect  observations  of  construction,  military 
and  amall  acale  field  explosive  events.  Where  liquefaction  has  been 
evidenced,  water  spouts,  sand  boils,  ground  subsidence  and  flow  of 
soils  have  followed  the  explosion  (Benson, 1983;  Blouin,1978;  Blouin  and 
Kim, 1983;  Carnes, 1976;  Charlie, 1978a, b;  Gilbert, 1976;  Melzer, 1978a) . 
It  is  also  believed  that  the  unusually  broad,  flat  crater  shapes 
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observed  ss  s  result  of  surface  and  near  surfaoe  detonations  in 
saturated  soils  stay  be  indicative  of  liquefaction  (Blouin,1978;  Blooin 
and  Kin, 19 83;  Carnes, 1976;  Gilbert, 1976) .  Explosive  loadings  pose  a 
twofold  potential  for  danage,  including  tbe  direct  effects  of  the 
explosion  as  well  as  the  additional  consequences  of  soil  liquefaction. 
Reviews  of  explosive-induced  liquefaction  phenomena  and  experience  are 
given  by  Charlie  et  al.  (1980,  1985),  Gilbert  (1976),  Marti  (1978)  and 
Rischbieter  (1977). 

Since  the  effects  of  liquefaction  from  explosive-induced  ground 
motions  are  not  restricted  to  surface  damage,  both  above  ground  and 
buried  structures  are  subject  to  potential  damage.  Therefore,  the 
possibility  of  explosive-induced  liquefaction  is  a  significant 
consideration  for  nearby  atructures.  The  resulting  damage  primarily 
involvea  bearing  capacity  failures,  including  excessive  settlements, 
upward  movements  and  rotations.  fhile  such  occurrences  may  not 
necessarily  damage  a  facility  structurally,  as  was  observed  in  the 
Niigata  earthquake,  they  may  render  it  unusable  from  a  practical  point 
of  view.  Therefore,  in  examining  the  liquefaction  potential  of  a  given 
site,  it  is  essential  to  maintain  a  proper  perspective  on  site  specific 
field  conditions  and  structural  considerations.  This  is  especially 
true  for  high  explosive  events  in  areas  where  the  groundwater  table  is 
within  the  depth  of  influence  of  explosively  generated  stresses. 

2.  Definition  of  Liquefaction 

The  term  liquefaction  has  been  generally  associated  with  the 
response  of  loose,  saturated  sands  to  loading  conditions  that  induce 
large  strains  resulting  in  flow  slides,  or  to  describe  a  state  of 
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stress  existing  where  the  effective  stress  spproaches  zero  or,  when  the 


redaction  in  soil  shear  strength  is  substantial  enough  that  the  soil 


behaves  like  a  liquid.  One  of  the  first  attempts  at  quantifying  this 


phenomenon  was  made  by  Casagrande  (1936).  Since  that  time  there  has 


been  increased  awareness  and  interest  in  liquefaction,  particularly  by 


researohers  investigating  seismic  soil  response.  Reviews  of  test 


procedures  for  dynamically  loading  soils  and  liquefaction  analyses  are 


given  by  Finn  (1972,1981)  and  Whitman  (1969). 


From  these  investigations,  a  number  of  definitions  and  terms 


associated  with  liquefaction  have  been  developed  (Committee  on  Soil 


Dynamics, 1978,  Seed, 1976).  Several  definitions  of  liquefaction  and 


related  observations  are  as  follows. 


Liquefaction:  'Denotes  a  condition  where  a  soil  will  undergo 

continued  deformation  at  a  low  residual  stress  or  with  no 
residual  resistance  due  to  the  build-up  of  high  porewater 
pressures  which  reduce  the  effective  confining  pressure  to  a 
very  low  value;  pore  pressure  build-up  leading  to  true 
liquefaction  of  this  type  may  be  due  to  either  static  or  cyclic 
stress  applications.'  (Seed, 1976). 


itlal  Liquefaction:  'Denotes  a  condition  where,  during  the 
course  of  cyclic  stress  application,  the  residual  porewater 
pressure  on  completion  of  any  full  stress  cycle  becomes  equal 
to  the  applied  confining  pressure;  the  development  of  initial 
liquefaction  has  no  implications  concerning  the  magnitude  of 
the  deformations  which  the  soil  might  subsequently  undergo; 
however  it  defines  a  condition  which  is  a  useful  basis  for 
assessing  various  possible  forms  of  subsequent  soil  behavior.' 
(Seed. 1976). 


allc  Liquefaction:  'Denotes  a  condition  in  which  cyclic  stress 
applications  develop  a  condition  of  initial  liquefaction  and 
subsequent  cyclic  stress  applications  cause  limited  strains  to 
develop  either  because  of  the  remaining  resistance  of  the  soil 
to  deformation  or  because  the  soil  dilates,  the  pore  pressure 
drops,  and  the  soil  stabilizes  under  the  applied  loads.' 
(Seed, 1976) . 
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Liquefaction:  'The  act  or  procesa  of  transforming  any  anbatance 
into  a  liquid.  In  cob*  a  ionic  a  a  soil,  the  transformation  ia 
froa  a  solid  atate  to  a  liquefied  state  aa  a  consequence  of 
increased  pore  pressure  and  reduced  effective 
stress. .. .Liquefaction  is  thus  defined  as  a  changing  of  states 
which  is  independent  of  initiating  disturbance  that  could  be 
static,  vibratory,  sea  wave,  or  shock  loading,  or  a  change  in 
groundwater  pressure.  The  definition  also  is  independent  of 
deformation  or  ground  failure  movements  that  might  follow  the 
transformation.  Liquefaction  always  produces  a  transient  loss 
of  shear  resistance  but  docs  not  always  produce  a  longer-term 
reduction  of  shear  strength.'  (Committee  on  Soil 
Dynamics. 1978) . 


Since  the  work  presented  herein  investigates  the  phenomenon  of 
liquefaction  as  a  result  of  impact  loadings  in  terms  of  the  porewater 
pressure  response  relative  to  the  initial  effective  stress  and  initial 
soil  density,  the  definition  by  Seed  (1978)  for  'initial  liquefaction' 
and  that  for  'liquefaction*  by  the  Committee  on  Soil  Dynamics  (1978) 
will  be  used  in  this  paper.  In  this  way  liquefaction  can  be  evaluated 
by  comparing  the  amount  of  porewater  pressure  increase  that  occurs 
after  the  attenuation  of  an  explosive-induced  or  shock-induced 
compressive  stress  wave.  The  porewater  pressure  increase  above  the 
hydrostatic  pressure.  remaining  after  the  dissipation  of  the 
compressive  stress  wave,  will  be  referred  to  as  the  'residual  excess 
porewater  pressure.' 

3.  Effective  Stress  and  Shear  Strenath 

The  behavior  of  soils  under  loads  is  governed  by  the  principle  of 
effective  stress  (Lambe  and  Whitman,19(9.  Terzaghi  and  Peck.1948). 
This  concept  was  first  presented  by  Terzaghi  (1943)  and  is  fundamental 
to  understanding  the  response  of  soils  under  various  loading 
conditions.  The  effective  stress,  o’  .  is  defined  as: 


where*  o  ie  the  totel  stress  end  u  is  the  porewater  pressure. 

The  shear  strength  of  a  soil.  t  »  proportional  to  the 
intergranular  pressure  and  the  internal  frictional  resistance  which  is 
expressed  by: 


o*  +  o'  Tan 


(Eq.  2.2) 


where,  c’  is  the  effective  cohesion  and  p'  is  the  effective  angle  of 
internal  friction.  For  cohesionless  soils  the  effective  cohesion,  c'  • 
is  xero  by  definition. 

The  condition  of  liquefaction  requires  a  complete  loss  of  shear 
strength  as  a  result  of  a  reduction  in  the  effective  stress  to  xero. 
The  effective  stress  is  xero  when  the  porewater  pressure  increases  in 
magnitude  equal  to  the  value  of  the  total  stress.  It  is  this  increase 
in  porewater  pressure  from  dynamic  loadings  which,  whether  from 
explosives,  earthquakes  or  shock  loadings,  can  induce  liquefaction  in  a 
saturated,  cohesionless  soil. 


B.  LlwfigUtm  Hffgfctalm 

i.  tolgiiiiziaflMifl  BsasaA  Matiam 

An  explosive  event  produces  high  amplitude,  high  frequency  ground 
motiona  of  relatively  short  duration.  The  intense  compressive  body 
waves  generated  by  an  explosion  propagate  outward  from  the  source  in  a 


cylindrieally  or  spherically  diverging  pattern,  depending  on  the 
geometry  of  the  source.  These  waves  attenuate  with  radial  distance,  R, 


froa  th«  diitubiBM.  For  cylindrical  wave  fronts  the  attenuation 
occurs  according  to  1/(R*^)  and  for  spherical  wave  fronts  according  to 

1/R  (Rinehart, 1975) .  Surface  waves  are  also  generated  and  propagate 
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froa  the  source  for  soae  distance,  attenuating  with  according  to  1/(R  ) 
(Richart  et  al.,1970).  Shear  waves  nay  also  be  produced  by  the 
reflection  of  coapression  waves  from  interfaces  such  as  the  ground 
surface,  groundwater  table,  soil  layers  or  underlying  bedrock.  Fron 
explosive  field  ground  notion  reoords,  it  is  evident  that  coapression, 
shear  and  surface  waveforns  exist  but  their  exact  interrelationship  in 
producing  liquefaction  ia  not  clear. 

Near  a  contained  (no  crater)  explosion,  the  soil  is  loaded  under 
confined,  uniaxial  conpreasion.  This  type  of  loading  has  been  observed 
in  strain  path  plots  of  field  data  fron  deeply  buried  tests  with 
spherically  shaped  charges  (Worknan  et  al.,  1981).  Radial  strains 
doninate  and  hoop  strains  are  virtually  nonexistant  until  after  the 
coaipreasive  stress  wave  has  sufficiently  attenuated.  When  tensile  hoop 
strains  are  developed  in  the  soil,  relative  grain  novenents  occur  and 
nay  disrupt  the  soil  fabric.  If  the  nagnitude  of  an  explosively 
generated  stress  wave  is  large  enough,  the  soil  will  fail  in  shear  as 
the  particles  separate  and  lose  grain-to-grain  contact 
(Melxer, 1978a, b) .  This  conplete  loss  of  shear  strength  is  consistent 
with  the  requirenent  for  liquefaction,  but  this  approach  does  not 
account  for  the  water  in  the  pores  or  the  snail  hoop  strains  at 
distances  away  fron  the  explosion. 

The  applied  conpressive  load  also  produces  a  direct  and  iamediate 
porewater  pressure  response  in  a  saturated  soil.  Correspondingly,  the 
soil  grains  and  soil  skeleton  are  also  stressed,  but  for  a  saturated 


•oil  the  pomittt  carries  a  large  part  of  the  load  increase.  Shock 
pulses  generated  by  explosions  typieally  result  in  nndrained  loading 
conditions.  particularly  near  the  source,  since  the  porewater  pressures 
cannot  diasipate  until  after  the  compressive  stress  wave  attenuates  and 
the  system  unloads.  During  unloading,  the  pressure  in  the  grain  matrix 
falls  abort  of  the  static  value  due  to  hysteresis  and  has  the  effect  of 
adding  to  the  porewater  pressure  which  reduces  the  effective  stress. 
The  onset  of  liquefaction  probably  occurs  during  the  unloading  phase  as 
the  systea  rebounds  (Schaeperaeier. 1978a) .  The  residual  excess 
porewater  pressure  decreases  the  effective  stress  reduces  the  soil's 
shear  strength.  Vith  tiae.  the  soil  gradually  regains  strength  as  the 
residual  exceas  porewater  pressure  dissipates  and  the  effective  stress 
is  restored. 

Seisaic  loading  of  the  soil  occurs  in  regions  where  the  compressive 
stress  waves  have  attenuated  and  the  ahear  waves  dominate.  In  these 
regions  liquefaction  aay  occur  froa  repeated  shearing  strain  reversals 
in  a  aanner  siailar  to  that  observed  for  earthquake-induced  ground 
aotiona.  For  surface  or  near  surface  explosions.  the  stress 
distributions  in  the  soil  becoae  more  complicated  by  the  presence  of 
air  overpressure  stresses  generated  by  the  exploaion  (Nelser. 1978b) . 
However*  in  soae  cases  the  near  surface  soils  may  undergo  a  one- 
diaensional  coapressive  loading  from  the  air  overpressure.  This 
condition  aay  be  analyzed  by  the  same  approach  suggested  above  for 
direct-induced  explosive  loadings. 


A  itinic  event  produce e  low  amplitude,  low  frequency  ground 
■otione  of  relatively  long  durations  and  nuaerout  cycles.  Shear  waves 
generated  during  an  earthquake  occur  as  a  result  of  stress 
redistributions  within  the  earth's  crust  that  release  energy 
proportional  to  the  fault  else  and  shear  aodulus.  Shear  waves  are 
generally  assuaed  to  propagate  upward  froa  the  underlying  bedrock 
through  the  soil  aass  in  an  oscillatory  aanner  involving  a  series  of 
stress  reversals.  Under  this  action,  the  soil  grains  begin  to 
rearrange  theaselves  relative  to  one  another  in  an  atteapt  to  reduce 
the  voluae  of  the  soil  skeleton  and  reach  a  condition  of  stable 
equilibrium  As  a  result  of  the  corresponding  reductions  in  fluid 
filled  pore  spaces,  the  porewater  pressure  increases  and  the  effective 
atress  decreases.  The  oscillatory  nature  of  the  earthquake  loading 
causes  the  porewater  pressure  to  gradually  increase.  Therefore,  the 
duration  and  nuaber  of  cycles  are  significant  factors  in  producing 
liquefaction  (Seed  and  Idriss.1971).  Since  the  porewater  is  not 
capable  of  aaintaining  any  shear,  the  shear  stresses  are  carried  by  the 
individual  grains  of  the  soil  skeleton.  When  liquefaction  occurs,  the 
intergranular  friction  is  teaporarily  lost,  and  the  soil  is  not  capable 
of  aaintaining  any  shear  strength.  With  tine,  the  residual  excess 
porewater  pressure  dissipates  through  fluid  drainage,  and  the  effective 
stress  and  shear  strength  are  gradually  regained. 

3.  Coanarison  of  Liquefaction  Mechanises 

Froa  the  nature  and  geoaetry  of  explosive  loadings,  both 


coapressive  and  shear  waves  are  generated.  The  coapressive  stress 
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waves  are  aost  significant  near  the  explosion  and  generate  intense 
coapressive  pulses.  Shear  stress  waves  doainate  at  distances  away  froa 
the  explosion  and  produce  oscillatory  notions.  Seisaic  ground  notions, 
siailar  to  those  produced  by  earthquake  loadings,  have  been  simulated 
in  the  field  by  controlled  explosive  tests  (Bruce  et  al.,  1979,1982). 
The  effect  of  these  stress  waves  in  producing  liquefaction  is  relative 
to  their  position  and  intensity  as  they  attenuate  away  froa  the  source. 

It  is  not  clear  how  the  unloading  and  shearing  actions  are  related 
near  the  explosion  and  it  is  difficult  to  precisely  define  the 
transition  between  then  froa  field  ground  notion  records.  However, 
these  actions  are  consistent  with  observations  of  relatively 
instantaneous  liquefaction  near  the  explosive  source.  At  distances 
away  froa  the  explosion,  seisaic  ground  notions  doainate  and 
liquefaction  aay  occur  in  a  aanner  siailar  to  that  froa  earthquake 
loadings  at  soae  tine  after  the  explosion  has  occurred.  The  resulting 
cyclic  shearing  action  gradually  increases  the  porewater  pressure  until 
the  effective  stress  goes  to  zero  and  the  soil  fails.  These 
observations  have  been  aade  by  a  number  of  investigators  (Benson, 1983; 
Charlie. 197 8b;  Gilbert, 1978;  Langley  et  al.,1972;  Keizer, 1978a) . 


C.  Pactora  Influenoina  Liquefaction 

There  are  a  significant  nuaiber  of  factors  that  have  been  observed 
to  influence  the  occurrence  of  liquefaction.  Most  of  the  experience 
with  exploaive  loadings  comes  indirectly  in  the  sense  that  soil 
densif ication  was  being  attempted  or  field  blasting  programs  were  being 
conducted  and  liquefaction  or  liquefaction  features  may  have  been 
noted.  The  general  observations  of  many  researchers  involved  in  both 


a 
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laboratory  and  field  experiments  investigating  compressional 
liquefaction  from  explosive  loadings  will  be  presented.  Some 
references  from  earthquake  researchers  will  also  be  given  where  similar 
observations  have  been  made. 

i.  foii-gtm.ar-t.tga. 

a.  Petree  of  Saturation 

Where  the  degree  of  saturation  is  less  than  100  percent* 
liquefaction  is  difficult  to  induce.  Even  small  amounts  of  air  can 
significantly  reduce  the  bulk  modulus  of  the  soil-water  mixture  and  the 
compressive  stress  wave  propagation  velocity.  Riohart  et  al.(1970) 
used  the  equations  for  the  compressibility  of  mixtures  and  the  wave 
propagation  velocity  in  mixtures  to  demonstrate  that  the  presence  of 
0.01  percent  air  in  a  aaturated  quarts  sand  will  reduce  the  bulk 
modulus  of  the  soil-water  system  by  a  factor  of  16.  aad  reduce  the  wave 
propagation  velooity  by  a  factor  of  four. 

While  liquefaction  of  materials  containing  air  bubbles  may  be 
difficult  to  induce  by  seismic  loadings,  it  may  be  possible  for 
liquefaction  to  occur  from  intense  shock  loadings  associated  with 
explosions.  This  can  occur  if  the  magnitude  of  the  applied  compressive 
stress  or  strain  is  large  enough  to  force  the  air  into  solution  and 
saturate  the  soil.  These  observations  have  been  made  by  a  number  of 
investigators  (Allen, 1975s  Allen  ct  si.. 1980;  Damitio, 1978b;  Florin  and 
Ivanov. 1961;  Perry, 1972;  Gilbert, 1976;  Ivanov, 1967;  Silver, 1981;  Studer 
aad  Kok,1980;  True ,1967, 1969;  Van  der  Kogel  et  al.,1981). 


b.  itltlill  PemitT 


Liquefaotion  susceptible  soils  beve  s  significant  potential 
for  volume  decrease.  Earthquake-induced  liquefaction  typically  occurs 
in  soils  baying  relative  densities  of  less  than  65  percent.  Loose 
soils  have  been  successfully  increased  in  relative  density  by  as  much 
as  40  percent  by  explosive  loadings  which  generally  results  in  a  aore 
stable  soil  aass.  Soae  dense  soils  nay  even  becone  loose  and  this  has 
been  aainly  attributed  to  the  dilative  behavior  of  dense  soils  which 
exhibit  a  volune  increase  accoapanied  by  a  reduction  in  porewater 
pressure.  These  observations  have  been  made  by  a  number  of 
investigators  (Allen, 1975;  Allen  et  al.,1980;  Daaitio. 1978b; 
D'Appolonia.1968;  Drake, 1978;  Finn, 1972;  Florin  and  Ivanov, 1961; 
Gilbert, 1976;  Ishihara  and  Vatanabe,1976;  Ivanov, 1967;  Klohn  et 
al.,1981;  Kok, 1978a;  Kuameneje  and  Eide,1961;  Kurzeme, 1971;  Mitchell 
and  Katti,1981;  Riachbieter.1977;  Seed  and  Idriss,1967;  Studer  and 
Hunziker,1977;  Studer  and  Kok, 1980;  Studer  and  Prater, 1977;  Yanaaura 
and  Koga.1974). 

c.  PtCTfftfrilitY 

The  duration  of  the  liquefied  state  depends  on  the 
permeability,  coapressibility,  stratum  thickness  and  drainage  path 
length.  The  rate  of  dissipation  of  the  residual  excess  porewater 
pressure  is  primarily  controlled  by  the  soil  permeability  which  is  a 
function  of  soil  type,  density,  grain  size,  packing,  cementation  and 
fluid  properties.  High  permeability  soils,  such  as  coarse  sands  and 
gravels,  tend  to  drain  rapidly  and  reduce  the  residual  excess  porewater 
pressure  quickly  enough  to  inhibit  liquefaction  and  restore  the 
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effective  stress.  Low  permeability  soils,  such  ss  fine  silts  end  low 
sensitivity  clays,  are  not  usually  prone  to  liquefaction.  Coarse  silts 
and  medium  to  fine  sands  have  been  observed  to  be  liquefaction 
susceptible  as  they  are  not  capable  of  relieving  the  residual  excess 
porewater  pressures  in  a  short  enough  period  of  tine.  Surface  features 
generally  believed  to  be  indicative  of  liquefaction  include  sand  cones, 
sand  boils,  saall  geysers  and  water  spouts  that  result  from  the 
dissipation  of  residual  excess  porewater  pressures.  These  features 
have  been  observed  at  various  tines  after  the  explosion  at  various 
distances  away  from  the  detonation  point,  depending  on  the  soil 
permeability  and  other  factors.  These  observations  have  been  nade  by  a 
number  of  investigators  (Daaitio, 1978a;  Das, 1983;  Florin  and 
Ivanov. 1961;  Gilbert, 1976;  Ivanov, 1967;  Ivanov  et  si., 1981;  Kok, 1978a; 
Kurxeae,1971;  Prater, 1977;  Prugh,1963;  Seed  and  Booker, 1976;  Studer  and 
Hunxiker.1977;  Studer  and  Kok, 1980;  Studer  and  Prater, 1977; 
True. 1967. 1969). 

d.  CflmiilUfrllitT 

The  transient  and  long-tern  porewater  pressure  responses  to 
applied  compressive  stress  loadings  are  dependent  upon  the  amplitude  of 
the  applied  stress,  the  compressibility  of  the  soil  skeleton,  the 
compressibility  of  the  pore  fluid  (wster  and/or  air)  and  drainage 
conditions.  The  presence  of  entrspped  air  increases  the 
compressibility  of  the  pore  fluid,  reducing  the  transient  and  long-term 
porewater  pressure  increases. 

The  tendency  for  volume  decrease  is  a  necessary  requirement  for 
liquefaction.  In  a  saturated  soil,  a  volume  decrease  causes  an 


increase  in  porewater  pressure  and  a  corresponding  decrease  in 


effective  stress.  As  such,  water  saturated  soils  having  a  low  relative 


density  should  have  a  high  liquefaction  potential.  The  potential  for  a 


soil  to  decrease  in  voluae  also  depends  on  the  stiffness  of  the  soil 


skeleton  which  is  a  function  of  both  density  and  effective  stress. 


Soils  having  low  effective  stress  conditions  should  have  a  high 


liquefaction  potential. 


The  presence  of  entrapped  air  increases  the  conpressibility  of  the 


systea,  reducing  the  potential  for  liquefaction  and  increasing  the 


attenuation  of  the  coapressive  stress  wave.  In  a  partially  saturated 


soil.  the  tendency  for  voluae  decrease  is  resisted  by  the 


coapressibility  of  the  soil  skeleton,  the  entrapped  air  and  the  water. 


whereas  in  a  saturated  soil  the  resistance  is  priaarily  due  to  the 


porewater.  The  pore  air  bubbles  tend  to  act  as  a  'cushion*  as  they 


absorb  the  coapressive  stress  during  deforaation.  These  observations 


have  been  aade  by  a  nuaber  of  investigators  (Florin  and  Ivanov, 1961; 


Gilbert,  1976;  Ishihara  and  Watanabe,1976;  Ivanov, 1967;  Jackson  et 


al.,1980;  Lyakhov, 1961;  Richart  et  al.,1970;  Studer  and  Kok,1980;  Van 


der  Kogel  et  al.,1981). 


In  general,  increasing  the  static  overburden  pressore 


increases  the  effective  stress  and  decreases  the  potential  for 


liquefaotion.  With  increasing  effective  stress  the  stiffness  of  the 


soil  skeleton  is  increased  and  the  porewater  pressure  aust  reach  a 


higher  value  before  liquefaction  can  occur.  These  observations  have 


been  aade  by  a  nuaber  of  investigators  (Allen, 1975;  Allen  et  si., 1980; 
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Daaitio, 1978a, b;  Florin  and  Ivanov, 1961;  Ivanov, 196 7;  Mitchell  and 
Katti,1981;  Oriard,1976;  Rischbieter,1977;  Seed, 197 9;  Seed  and 
Idrisa,1971i  Studer.1978;  Stnder  and  Hunziker,1977;  Studer  and 
Prater, 1977). 

f.  Coheiion 

Observations  from  field  blasting  experience  indicate  that  the 
presence  of  cohesion  in  soils  tends  reduce  the  liquefaction  potential 
(Daaiitio, 1978b;  Gilbert, 1978;  Rischbieter ,1977) .  It  has  generally  been 
observed  that  granular  soils  with  little  or  no  cohesion  are  the  most 
susceptible  to  liquefaction  for  seismic  loadings.  D'Appolonia  (1968) 
suggests  that  less  than  10  percent  silt  or  clay  size  particles  is 
sufficient  for  a  granular  soil  to  be  considered  cohesionless.  Seed 
(1981)  indicates  that  soils  containing  more  than  10  percent  clays  are 
generally  not  sensitive  to  liquefaction  from  earthquake  loadings. 

g.  fllidltlPB 

It  has  been  recognized  by  most  researchers  that  for  explosive 
and  seisaie  loadings,  unifomly  graded  soils,  having  a  narrow  band  of 
grain  sizes,  tend  to  be  nore  liquefaction  susceptible  than  soils  having 
a  variety  of  particle  sizes.  The  distribution  of  particle  sizes  in  a 
well  graded  soil  are  effective  in  filling  pore  spaces  with  soil  grains 
which  increases  the  particle-to-particle  oontact  area  and  shearing 
resistance.  These  observations  have  been  siade  by  a  number  of 


investigators  (Damitio, 1978b;  D'Appolonia,  1968;  Klohn  et  al.,1981; 
Kuaaeaeje  and  Eide,1961;  Rischbieter, 1977;  Seed  and  Lundgren,1954; 


Studer, 1978;  Studer  and  Hunziker,1977;  Studer  and  Kok, 1980;  Studer  and 


Prater, 1977;  True, 1967, 1969) . 

Studer  and  Kok  (1980)  have  noted  substantial  porewater  pressure 
increases  in  coarse  granular  soil  under  shock  loadings  in  the 
laboratory.  Florin  and  Ivanov  (1961)  have  demonstrated  in  laboratory 
experiments  that  all  sufficiently  loose,  cohesionless  soils  of  any 
grain  size  can  be  liquefied  under  explosive  loadings. 

h.  Particle  Size  andShane 

Fine  sands  have  been  found  to  be  generally  the  most 
susceptible  to  liquefaction  as  compared  to  other  grain  sizes  for  both 
explosive  and  seismic  loadings  (Damitio, 1978b;  D'Appolonia,  1968; 
Gilbert, 1976;  Klohn  et  al.,1981;  Knmmeneje  and  Eide,1961;  Seed  and 
Lundgren,19S4;  Studer, 1978;  Studer  and  Hunziker,1977;  Studer  and 
Prater, 1977;  Yang. 1973). 

Studer  and  Kok  (1980)  have  noted  substantial  porewater  pressure 
increases  in  coarse  granular  soil  under  dynamic  loadings  in  the 
laboratory.  Florin  and  Ivanov  (1961)  have  demonstrated  in  laboratory 
experiments  that  all  sufficiently  loose,  cohesionless  soils  of  any 
grain  size  can  be  liquefied  under  explosive  loadings. 

There  is  relatively  little  information  available  on  particle  shape. 
True  (1967,1969)  observed  that  a  graded,  angular  sand  showed  more 
evidence  of  liquefaction  than  a  uniform,  rounded  sand  under  shock 
loadings  in  the  laboratory.  Liquefaction  of  crushed  quartz  sand  was 
observed  by  Studer  (1978)  and  Studer  and  Prater  (1977)  for  explosive 
and  shock  loadings.  Klohn  et  al.  (1981)  were  able  to  liquefy 
subangular  tailings  by  explosions.  Ivanov  et  al.  (1981)  observed  that 
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angular  to  rabtaplir  cohesionless  sands  had  a  store  stable  structure 
under  dynaaic  loadings  than  did  rounded  cohesionless  sands.  Based  on 
the  United  data  available,  it  appears  that  it  is  possible  to  liquefy 
both  angular  and  rounded  cohesionless  sands. 

i.  BUlUm 

Dilatency  is  believed  to  play  an  important  role  in  the  early 
tine  response  of  explosively  loaded  soils  having  high  initial  densities 
(Drake, 1978;  Gilbert, 1976;  Lyakhov, 1961;  Seed, 1979;  Seed  and 
Lundgren,1954;  Studer  and  Kok,1980).  In  nediun  dense  and  dense  sands 
the  dilatent  nature  of  these  states  tends  to  inhibit  the  onset  of 

liquefaction.  As  shearing  ocours  in  the  soil  nass,  the  grains  slide 
over  and  past  one  another  causing  a  volune  expansion  of  the  soil 

skeleton  (Laabe  and  Vhitaian,1969) .  The  rearrangeaent  of  the  soil 

skeleton  results  causes  an  increase  in  void  sixes  that  dissipates  the 
residual  excess  porevater  pressures  and  increases  the  effective  stress. 
In  loose  soils  there  is  a  tendency  for  a  voluaie  decrease  to  occur  in 
the  soil  skeleton  as  the  soil  grains  try  to  rearrange  theaselves  into  a 
denser  packing  arrangeaent.  This  rearrangeaent  reduces  the  void  sizes 
which  causes  the  residual  excess  porevater  pressure  to  increase  and  the 
effective  stress  to  decrease.  Casagrande  (1936)  described  this 
behavior  and  atteapted  to  quantify  it  by  noting  that  there  is  a 

'critical  density'  at  vhich  no  voluae  change  will  occur. 

j.  Local  Geoloav 

The  presence  and  location  of  the  groundwater  table,  along  with 
geologic  interfaces  including  underlying  bedrock,  often  produce 


secondary  oscillatory  shear,  tensile,  compressive  and  surface  waves 
that  may  influence  and  enhance  the  possibility  of  liquefaction  under 
explosive  loadings.  Stiff,  dense  soils  and  rock  transmit  high 
frequency  and  high  amplitude  ground  motions.  The  magnitude  of  these 
ground  motions  depends  on  the  proximity  of  the  soil  and  rock  to  the 
explosive  detonation  point  (Drake  and  Ingram, 1981;  Ivanov, 1967;  Oriard, 
1976;  Sanders, 1982) . 

k.  gyjJm&iit 

When  material  damping  is  high,  the  attenuation  of  explosively 
generated  stress  wave  energy  is  greater  and  the  extent  of  the  liquefied 
zone  is  reduced.  Damping  is  a  function  of  frequency,  compressibility, 
strain  amplitude  and  material  properties.  The  denser  and  stiffer  the 
soil  or  rook  the  higher  the  frequency  and  amplitude  of  the  transmitted 
explosive-induced  ground  motions  (Oriard,1976;  Richart  et  al.,1970; 
Rinehart,  1975).  Dynamic  tests  of  dry  granular  soils  show  that  they 
exhibit  an  'S'-shaped  stress-strain  curve  and  a  locking  phenomenon.  It 
is  believed  that  this  may  have  a  significant  effect  on  the  rise  time  of 
explosive-induced  stresses,  involving  a  strain  rate  soil  strength 
dependency  (Jackson  et  al.,1980;  Whitman, 1957;  Whitman  and  Healy,1962; 
Whitman  et  al.,1964). 


2 .  Bxploaive  Loading  Parameters 
a.  Lowfl  lalwaiJJa 

The  potential  for  liquefaction  is  strongly  influenced  by  the 
magnitude  of  the  applied  stress  and  resulting  strains.  Larger 
intensity  charges  are  necessary  to  induce  liquefaction  with  increasing 
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depth  and  initial  density.  A  series  of  consecutive  explosions  are 
usually  aore  effective  due  to  the  superposition  of  applied  stresses. 
These  observations  have  been  nade  by  a  number  of  investigators 
(Daaitio, 1978a;  Florin  and  Ivanov. 1961;  Ball, 1962;  Ivanov, 1967;  Klohn 
et  al.,1981;  Mitchell  and  Katti,1981;  Studer  and  Prater, 1977;  Yamamura 
and  Koga,1974;  Yang, 1973). 

b.  Depth  of  Burial  of  Charae 

A  contained  explosive  charge  is  aore  effective  in  producing 
liquefaction  than  a  surface  charge  of  comparable  weight  and  density 
since  aore  of  the  explosive  energy  is  transaitted  to  the  soil. 
Explosives  are  aost  effective  in  producing  liquefaction  when  placed  at 
a  depth  of  about  two  thirds  the  stratus  thickness  of  the  potentially 
liquefiable  soil  deposit  (Daaitio, 1978a;  Hall, 1962;  Mitchell, 1970; 
Mitchell  and  Katti.1981;  Prug h,1963) . 

Under  given  circumstances  (i.e.,  loose,  saturated,  cohesionless 
soil  within  the  influence  of  a  sufficiently  intense  explosive  loading), 
the  energy  released  froa  a  surface  or  near  surface  explosion  of 
sufficient  intensity  can  produoe  liquefaction  (Charlie, 197 8b; 
Benson, 1983;  Blouin,1978;  Blouin  and  Kin, 1983;  Daaitio, 1978a, b; 
Gilbert, 1976;  Ivanov, 1967;  Yang, 1973).  It  is  also  feasible  that  a 
surface  and  a  buried  charge  can  be  scaled  accordingly  to  produce  the 
saae  depth  of  liquefaction.  Using  the  eapirical  relationships 
presented  in  the  later  part  of  this  chapter,  the  required  surface 
charge  is  estiaated  at  about  10  tines  that  required  for  a  buried  charge 
to  obtain  the  saae  depth  of  liquefied  soil. 


A«  the  density  of  e  charge  or  series  of  charges  increases,  the 
ability  to  induce  liquefaction  also  increases  since  a  load  of  greater 
intenaity  will  be  prodnced.  Experience  has  shown  that  the  most 
efficient  charge  shape  is  a  sphere.  However,  since  the  weight  mast 
increase  with  depth  of  bnrial  to  be  effective,  cylindrical  charges 
having  a  height  of  abont  three  to  four  times  the  diameter  are  generally 
used  (Damitio, 1978a;  Ivanov, 1967;  Oriard,197 6). 

d.  Charae  Weiaht 

In  general,  ground  shock  amplitude  increases,  frequency 
decreases  and  the  extent  of  the  liquefied  xone  increases  with 
increasing  charge  weight.  Repeated  loadings  with  small,  properly 
sequenced  chargee,  are  the  moat  effective  in  producing  liquefaction. 
However,  single  large  charges  of  sufficient  intensity  are  also  capable 
producing  liquefaction  (Charlie, 1978b;  Benson, 1983;  Blouin.1978;  Blouin 
and  Kim,1983;  Damitio, 1978a, b;  Gilbert, 1976;  Hall, 1962,  Ivanov, 1967; 
Lyakhov, 1961;  Oriard,1976;  Tang, 1973). 

••  PpIiy  XtoBtia.  fluu  IAllPxai-hajLMi«hax  af  £x&la± 

Experience  has  shown  that  properly  sequenced,  successive 
explosions  are  more  effective  in  producing  liquefaction  than  a  single, 
large  explosion.  The  time  delayed  loadings  take  advantage  of  the 
superposition  of  applied  stresses  and  the  effects  of  repeated  loadings 
(Damitio,1978a;  Hall, 1962;  Ivanov, 1967;  Kummeneje  and  Eide,1961; 


Orisrd.1976) 


For  explosive  loadings,  the  duration  of  the  applied  coapressive 
stress  wave  the  superposition  of  stresses  f  roa  consecutive  delayed 
detonations  ean  induce  liquefaction.  If  the  porevater  pressures  are 
allowed  to  dissipate  between  detonations,  the  soil  can  stabilise  and 
aay  becoae  aore  resistant  to  further  liquefaction.  These  observations 
have  been  aade  by  a  nuaber  of  investigators  (Daaitio, 1978a, b;  Florin 
and  Ivanov, 1961;  Ivanov, 1967;  Kuaaeneje  and  Eide,1961;  Oriard,1976; 
True, 1967. 1969). 

The  geoaetrical  attenuation  of  stress  waves  propagated  during  an 
explosive  event  can  be  altered  by  the  proper  selection  of  delay  tiaes 
and  charge  patterns  to  produce  stress  distributions  that  can  induce 
liquefaction.  Superiaposed  shock  waves  f roa  siaultaneous  explosions 
can  increase  the  extent  of  the  liquefied  xone.  However,  the 
appropriate  sequencing  of  delay  tiaes  presents  a  coaplicated  problea. 
If  the  delay  tiaes  between  successive  charge  detonations  is  too  short, 
then  incoaplete  liquefaction  occurs.  If  the  delay  tiaes  are  too  long, 
then  the  soil  is  allowed  to  drain  and  stabilize  before  being  fully 
liquefied  and  aay  in  fact  becoae  aore  resistive  to  liquefaction.  These 
observations  have  been  aade  by  a  nuaber  of  investigators 
(Daaitio, 1978a, b;  Hall, 1962;  Ivanov, 1967;  Klohn  et  si., 1981;  Kuaaeneje 
and  Bide, 1961;  Oriard,1976;  Mitchell, 1981;  Prugh.1963). 

The  nuaber  of  cycles  and  the  period  of  application  are  significant 
factors  in  producing  liquefaction  for  seisaic  loading.  For  exaaple,  if 
the  duration  of  the  1964  Alaskan  earthquake  had  been  shorter, 
liquefaction  would  not  have  occurred  over  such  an  extensive  area  (Seed 
and  Idriss,1971) .  The  superposition  of  successive  explosive 
detonations  also  produces  a  series  of  stress  wave  'cycles'  that  can 


cause  liquefaction  if  the  residual  excess  pore water  pressure  is  not 
allowed  to  dissipate  between  charge  delays  (Danitio  1978a). 


f.  Qe one trie  Attenuation 

The  attenuation  of  the  applied  oonpressive  stress  wave  depends 
on  aany  factors  including  the  stiffness  and  density  of  the  soil  or 
rock,  soil  saturation,  soil  compressibility,  charge  geometry,  load 
intensity  and  duration  of  the  applied  load.  The  geometric  spreading  of 
stress  waves  accounts  for  a  major  portion  of  the  attenuation  with 
distance.  Cylindrically  shaped  charges  generate  stress  waves  that 
propagate  from  the  source  in  a  cylindrically  diverging  pattern. 
Spherically  shaped  charges  generate  stress  waves  that  propagate  from 
the  source  in  a  spherically  diverging  pattern.  Both  oharge  geometries 
attenuate  the  stress  wave  in  an  ever  expanding  radial  manner.  The 
amount  of  attenuation  generally  tends  to  inorease  with  depth  and 
distance  away  from  the  detonation  point.  These  observations  have  been 
made  by  a  number  of  investigators  (Drake  and  Ingram, 1981;  Klohn  et 
al. ,1981;  Lyakhov. 1961;  Nelzer ,1978b;  Oriard.1976;  Rinehart, 1975; 
Rischbieter  et  al.,1977;  Sanders, 1982;  True, 1967, 1969) . 


3.  Swiff TT 

As  can  be  seen  by  the  number  of  parsmeters  that  influence  the 
occurrence  of  liquefaction,  a  complete  characterization  of  a  given  site 
becomes  extremely  difficult  and  complex.  Evaluating  each  parameter  is 
a  cumbersome  task  with  much  uncertainty.  There  is  some  control  that 
can  be  exercised  in  connection  with  the  selection,  location  and 
detonation  of  the  explosive  charge,  but  this  is  limited  in  terms  of 


sit*  response.  Under  such  circumstances,  it  is  necessary  to  treat 
these  factors  in  a  general  manner  with  a  number  of  simplifying 
assumptions  so  tbat  this  information  can  be  applied  to  a  variety  of 
sites.  Obviously  this  limits  the  applicability  of  general  trends 
observed  to  different  soil  types.  However,  for  any  particular  case,  it 
should  be  emphasised  that  site  speoific  conditions  must  be  thoroughly 
investigated  and  incorporated  into  a  comprehensive  analysis  and 
evaluation. 


The  detonation  of  high  yield  ezploaives  by  the  military  has  been 
done  primarily  to  investigate  structural  response  and  crater  formation 
processes  (Crawford  et  al.,  1974).  Considerable  attention  has  also 
been  given  to  evaluating  the  dynamic  in  situ  properties  of  various 
geologic  materials.  However,  little  consideration  has  been  given  to 
evaluating  the  liquefaction  phenomenon  for  these  events.  Available 
data  on  porewater  pressure  response  are  limited  at  best  and  are  often 
sketchy  and  incomplete.  Indirect  observations  of  liquefaction  have 
been  made  at  some  sites  snd  do  provide  some  useful  information.  At 
sites  having  cohesionless  soils  where  the  groundwster  tsble  is  nesr  the 
surface,  the  craters  are  typically  broad  and  flat  in  shape.  In  areas 
where  the  groundwater  table  is  not  within  the  depth  of  influence  of  the 
explosion,  craters  generally  assume  a  deep  bowl  shaped  form.  Such 
observations  have  only  recently  caused  a  re-evaluation  of  past  events 
where  unusual  cratering  phenomena  have  been  noted  and  it  is  believed 
this  could  possibly  be  explained  in  terms  of  liquefsetion. 
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At  the  Pacific  Proving  Grounds,  possible  indicators  of  explosive- 


indnccd  liquefaction  include  broad,  flat  craters,  sand  cones,  sink 


holes  and  other  collapse  features  beyond  the  crater  ria  (Blouin,1978; 


Blouin  and  Kim,1983;  Couoh  et  al.,1975;  Ristvet  et  el., 1978).  The 


SNOWBALL  event  also  had  a  broad,  flat  crater  with  water,  sand  and  silt 


flowing  into  the  area  after  the  explosion  (Gilbert, 1976;  Jones, 1976; 


Roddy, 1976a, b) .  Similar  trends  were  observed  at  the  PRAIRIE  FLAT 


event,  including  sand  boils,  water  spouts  and  ground  subsidence 


(Melxer, 1978a;  Roddy, 1976b) .  Piexoaeters  installed  at  the  DIAL  PACK 


event  test  site  indicated  that  residual  excess  porewater  pressures 


existed  for  soae  tine  after  the  explosives  were  detonated 


(Charlie, 1978a,  Langley  et  al.,1972).  The  PRE-DICE  THROW  II  event 


exhibited  siailar  behavior.  Springs  flowed  for  several  hours  after  the 


explosion  and  the  resulting  crater  depth  was  about  half  that  predicted 


(Helxer,1978a) .  The  DISTANT  RUNNER  II  field  test  foraed  a  large,  deep. 


bowl  shaped  crater  initially,  with  water  later  flowing  in  and  nearly 


filling  it.  The  DISTANT  RUNNER  III  event  produced  a  wide,  flat  crater 


that  filled  with  water  and  sand  immediately  after  the  explosion.  The 


DISTANT  RUNNER  events  occurred  about  3S0  meters  apart  in  similar 


geologic  materials,  but  had  markedly  different  cratering 


characteristics.  The  differences  in  the  resulting  crater  formations  is 


believed  to  be  due  to  variations  in  the  soil  profiles  and  liquefaction 


effeots  at  these  two  sites  (Benson, 1983 ) . 


All  these  events  took  place  in  loose,  saturated  soils  where  the 


groundwater  table  was  shallow  relative  to  the  depth  of  influence  of  the 


explosion.  It  should  be  noted  that  except  for  the  DIAL  PACK  event. 


there  is  virtually  no  data  available  from  these  sites  for  any  direct 


HiiUMenti  of  porewater  prostare  response.  Considering  thet  soil  end 
groundwater  eonditions  for  the  oeoarrenee  of  liquefaction  exist  in  aany 
areas  of  the  world,  particularly  when  snbjeoted  to  explosive  loadings, 
this  lack  of  data  is  unfortunate. 

2.  Conatruotion  Blastina 

Much  of  the  field  experience  in  explosive-induced  liquefaction 
has  coae  f row  the  practical  application  of  explosives  to  engineering 
probleas.  As  early  as  1932.  explosives  were  used  to  coapaot  fill  in 
the  United  States  (Engineering  News  Record. 1932) .  Since  that  tiae, 
there  have  been  aany  reports  of  field  blasting  prograas  involving  the 
densif ication  of  loose  sand  deposits,  stabilisation  of  subaerged  slopes 
and  site  preparation  (Abelev  and  Askalonov.1957;  Ivanov, 1967; 
Hall, 1962;  Kunmneje  and  Eide.1961;  Klohn  et  al.,1981;  Kok, 1978c;  Long 
et  al.,1981;  Lyaan.1942;  Mitchell  and  Katti.1981;  Mueller, 1971; 
Obemeyer.1980;  Prakash  and  Gupta, 1970;  Prugh.1963;  Quieroz  et 
al.,1967;  Sanders, 1982;  Solyaar,1984;  Solyaar  et  al.,1984; 
Terxaghi,1956;  Wild  and  Haslaa,1962) .  During  these  events,  care  is 
often  taken  to  aonitor  the  porewater  pressure  response  and  settleaents 
and  to  tailor  the  field  blasting  prograa  to  ainiaize  the  potential  for 
liquefaction.  These  operations  have  provided  soae  data  and  practical 
observations  of  soil  and  porewater  pressure  response  under  field 
explosive  loadings. 

3.  8—11  Se«ie  Field  Exneriaents 

Soae  work  has  been  done  with  field  experiaents  specifically 
designed  to  investigate  soil  liquefaction  and  coapaction.  Russian  and 


European  researchers  have  perforated  tests  in  saturated,  sandy  soils  to 
investigate  the  effect  of  charge  weight,  location,  patterns  and  delay 
times  on  liquefaction  (Ivanov, 1967;  Kok, 1977 , 1978a, b;  Lyakhov, 1961; 
Puchkov, 19(2;  Rischbieter,1978;  Schaeperaeier, 1978b, o;  Studer  et 
al.,1978;  Trense ,1977) .  These  investigations  have  led  to  the 
developaent  of  eapirical  equations  that  can  be  used  to  estiaste  the 
extent  of  liquefaction  at  a  site  in  teras  of  explosive  charge 
paraaeters.  Siailar  investigations  have  been  conducted  by  researchers 
in  Japan  (laaaaura  and  Koga.1974)  and  India  (Arya  et  al.,1978;  Prakash 
and  Gupta,1970).  The  U.S.  Army  Corps  of  Engineers  Waterways  Experiment 
Station  has  conducted  field  explosive  tests  at  several  locations  to 
investigate  the  influence  of  charge  location  relative  to  the 
groundwater  table  (Drake, 1978;  Carnes,  1976;  Perry, 1972).  Evidence  of 
liquefaction  was  not  observed  where  dense,  coarse  or  very  fine  sands 
existed  or  where  conditions  of  saturation  were  not  evident.  It  has 
been  noted  that  where  broad,  shallow  craters  were  produced  in  these 
small  scale  tests,  loose,  saturated  soils  were  within  the  depth  of 
influence  of  the  explosion.  Conversely,  where  deep  bowl  shaped  craters 
were  formed  this  was  not  the  case.  These  observations  are  consistent 
with  experience  gained  in  military  high  explosive  events. 


E.  Laboratory  Experience 

1.  Small  Scale  Explosive  Tests 

A  limited  amount  of  tests  have  been  conducted  using  small  scale 
explosive  events  in  the  laboratory.  The  Bussians  have  investigated  the 
effects  of  surface,  buried  and  underwater  explosions  in  loose, 
saturated,  sandy  soils.  These  soils  were  tested  in  metallic  tanks  and 
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reinforced  count*  trays  (Ivsnov,1967) .  Soil  stress  and  porewater 
pressure  measurements  were  asde  end  the  influence  of  sncoessive 
detonations  and  delays  was  exaained.  Bnried  charges  and  appropriately 
sequenced  delays  were  found  to  be  the  most  effective  in  producing 
liquefaction.  This  is  consistent  with  results  obtained  in  saall  scale 
field  tests.  Laboratory  tests  have  also  been  conducted  to  evaluate  the 
effect  of  charge  placement  on  cratering  in  loose,  saturated,  sandy 
soils  (Piekutowski.l97<;  Vesie  et  si., 1967). 

Attempts  have  been  aade  to  use  the  centrifuge  to  aodel  field 
explosion  conditions  (Al-Hussani,1976;  Fragaszy  et  S1..1983; 
Janes, 1978;  Schmidt. 1976;  Schmidt  et  al.,1981).  This  laboratory 
technique  involves  the  subscale  modeling  of  prototype  stress  levels 
where  small  explosive  charges  are  detonated  under  the  influence  of  an 
acceleration  field.  Results  have  indicated  that  it  is  possible  to 
induce  liquefaction  by  this  method.  Some  measurements  have  been  made 
of  soil  and  porewater  pressure  response.  Cratering  observations  have 
generally  been  consistent  with  field  experience.  Fragaszy  et  al. 
(1983)  have  presented  their  findings  in  terms  of  a  porewater  pressure 
ratio  (the  ratio  of  the  residual  excess  porewater  pressure  to  the 
initial  effective  stress)  as  a  function  of  charge  weight  and  this 
relationship  is  presented  in  Figure  2.1. 

Experiments  have  been  performed  using  air  shock  loadings  to 
develop  compressive  stresses  on  the  order  of  0.20  MPa  with  millisecond 
rise  times.  A  direct  stress  loading  (Van  der  Kogel  et  al.,1981)  and  an 
air  overpressure  loading  that  passes  over  the  soil  have  been  used 
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(Studer,1978;  Studer  and  Prater, 1977;  Quaak.1978) .  The  direct  loading 
teats  were  done  on  dry  and  partially  saturated  sands.  Test  results 
indicated  that  rise  tines  for  porewater  pressure  response  decreased 
with  increasing  degrees  of  saturation.  The  air  overpressure  loadings 
were  done  on  loose,  saturated  sands.  Liquefaction  was  observed  as  a 
result  of  residual  excess  porewater  pressures  reaiaining  after  the 
stress  wave  was  passed  over  the  samples.  Perry  (1972)  investigated  the 
effects  of  density  and  aoisture  content  on  the  liquefaction  potential 
of  loose  sands  using  an  explosive  overpressure  simulation  device. 
Liquefaction  was  evaluated  by  monitoring  the  relative  movements  of  a 
number  of  variable  density  inclusions  in  the  sand  after  loading.  In  a 
few  oases  liquefaction  was  observed  at  low  relative  densities  and  high 
moisture  contents. 

Impact  loading  tests  of  saturated  and  nearly  saturated  sands  have 
also  been  done.  Tanimoto  (1987)  used  a  swinging  pendulum  to  impact  the 
side  of  a  vertical  column  of  loose,  saturated  sand.  Measurements  were 
made  of  porewater  pressure  increases  and  settlements  caused  by 
liquefaction.  Greater  settlement  occurred  with  greater  increases  in 
porewater  pressure.  Kok  (1977)  also  used  a  pendulum  to  apply  a 
horixontal  impact  to  a  vertical  column  in  order  to  observe  the  effect 
of  applied  stress  intensity  on  the  gradient  across  the  soil.  The  more 
intense  the  applied  load,  the  larger  the  gradient  across  the  sample. 
Ruygrok  and  Van  der  Kogel  (1980)  used  a  free  falling  weight  system  to 
simulate  the  loading  applied  by  ocean  waves  on  covered,  loose,  nearly 
saturated  sand  and  were  able  to  obtain  liquefaction  in  a  number  of 
tests.  Allen  (1975)  and  Allen  et  si. (1980)  investigated  the  effects  of 
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saturation  and  relative  density  on  stress  wave  propagation  velocity, 
observing  higher  velocities  for  dense,  saturated  soils. 

Dynamic  raa  loading  systems  have  been  developed  at  the  Vatervays 
Experiaent  Station  (Jackson, 1978;  Jackson  et  si., 1980;  Seaaan,1983) . 
These  devices  are  capable  of  generating  shock  loads  of  up  to  400  MPa 
with  Billisecond  rise  tines  on  uniaxial  and  triaxial  test  specimens. 
Dynaaic  ooapressibility,  vertical  stress  and  deformation  can  be 
measured.  Based  on  an  analysis  of  test  data,  it  has  been  shown  that 
strain  rate  effects  become  significant  for  submillisecond  rise  times  in 
dry  and  partially  saturated  soils.  Plans  are  being  made  at  this  time 
to  incorporate  a  dynamic  porewater  pressure  measurement  system. 

3.  Quasi-Statio  Loading 

Experiments  using  quasi-static,  isotropic  compressive  loadings 
have  been  performed  to  investigate  the  liquefaction  susceptability  of 
saturated  coral  sands  and  saturated  quarts  sands  (Fragaszy  and 
Voss, 1981;  Fragaszy  et  al.,1983).  The  coral  sand  was  obtained  from  the 
Eniwetok  Atoll  in  the  south  Pacific  which  was  the  site  of  nuclear 
explosive  detonations  in  the  1950's  (Couch  et  al.,  1975;  Ristvet  et 
al.,  1978).  A  complete  loading  cyole  was  applied  in  one  to  two  minutes 
and  applied  stresses  were  on  the  order  of  35  MPa.  Both  drained  and 
undrained  behavior  were  exaained  and  some  evidence  of  liquefaction  was 
observed  in  a  few  tests  during  unloading.  Results  indicate  that  grain 
fracturing  nay  be  a  significant  factor  in  liquefying  coral  sands  and 
that  it  is  possible  to  induce  liquefaction  by  applying  an  isotropic 
coapressive  loading.  Although  porewater  pressure  increases  were 


observed  in  the  quartz  sand,  liquefaction  did  not  occur. 


The  empirical  relationships  that  are  available  have  been  developed 

# 

from  field  experience  using  buried  charges  in  loose,  saturated, 
granular  soils.  They  are  useful  approximations  that  can  be  applied  in 
a  general  manner  to  a  variety  of  soils  in  field  blasting  programs  where 
buried  charges  are  used.  However,  they  should  be  used  with  appropriate 
judgement  and  consideration  for  particular  site  specific  conditions. 

1.  Explosive  Charge 

Baaed  on  field  experience  and  observations,  several  empirical 
relationships  have  been  developed  to  estimate  the  extent  of  the 
liquefied  xone  and  the  radius  of  influence  of  contained  charges.  These 
equations  were  obtained  using  buried,  contained  explosions  unless 
otherwise  indicated. 

The  required  charge  weight,  W,  in  kilograms  per  delay  to  insure  a 
contained  explosion,  is  given  by  (Ivanov, 1967;  Damitio, 1978b) : 

1  -  (0.55)  d3  (Eq‘  2,3> 


where,  d  is  the  buried  charge  depth  in  meters. 

The  depth  of  the  liquefied  soil,  1^  ,  in  meters  for  the  charge 
weight,  V,  in  kilograms  per  delay  given  in  Equation  2.3,  is  given  by 
(Ivanov,  1967;  Damitio, 1978b) : 

(Eq.  2.4) 

where,  d  corresponds  to  that  used  in  Equation  2.3. 


The  optimum  placement  depth,  h  ,  in  meters  for  a  charge  weight, 
V,  in  kilograms  per  delay  to  obtain  the  largest  radios  of  liquefaction 
at  the  ground  aurface,  is  given  by  (Lyakhov, 1961) : 


(2.5)  W 


1/3 


(Eq.  2.5). 


The  maximum  radius  of  the  liquefied  region,  R  ,  in  meters,  for 
the  optimum  placement  depth,  hQ  ,  in  meters  from  Equation  2.5,  is  given 


by  (Lyakhov, 1961) 


R 


max 


(b)  W 


1/3 


(Eq.  2.6) 


where,  V  is  the  charge  weight  in  kilogram)1;  per  delay  and  b  is  a 
coefficient  based  on  the  charge  weight  and  scaled  depth  of  burst. 


d/1 


1/3 


,  given  in  Figure  2.2. 


The  maximum  radius  of  the  liquefied  region,  R _  ,  in  meters  for 

max 


the  charge  placement  depth  used  in  Equation  2.3,  is  given  by 
(Ivanov, 1967;  Damitio, 1978b) : 


R 


max 


k3  W 


1/3 


(Eq.  2.7) 


where,  k^  is  an  empirical  factor  based  on  relative  density  and  grain 
size,  having  values  given  in  Table  2.1. 

% 

The  effective  radius  of  influence  of  a  buried  charge,  R 
meters,  is  given  by  (Ivanov, 1967;  Damitio, 1978b) : 


eff 


in 
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Table  2.1  Empirical  Sealing  Faetora.  k.;.  and  k  ,  for  Use 
for  Dae  in  Equations  2.7  and *2.8 
(Ivanov, 1967;  Daaitio, 1978b) 


Type  of  Soil 

Relative  Density 

k3 

k4 

Fine  Sand 

0.0  -  0.2 

25  -  15 

5-4 

Fine  Sand 

0.3  -  0.4 

9-8 

3 

Fine  Sand 

>  0.4 

<  7 

<  2.5 

Medium  Sand 

0.3  -  0.4 

8-7 

3  -  2.5 

Medium  Sand 

>  0.4 

<  6 

<  2.5 

•is*,  having  values  given  in  Table  2.1.  In-plane  charge  spacings 
should  be  limited  to  two  tines  the  effective  radius  of  influence  for 


naxiaun  effectiveness. 

An  estimate  of  the  maximum  depth  of  liquefaction.  h(  ,  in  meters. 
«as  determined  from  tests  using  surface  detonations  of  explosives  and 
is  expressed  by  Equation  2.9.  Due  to  the  lover  ground  shock  energy 
transmission  from  surface  explosions,  the  extent  of  the  liquefied  zone 
is  less  than  for  buried  charges  (Ivanov. 1967;  Daaitio. 1978b) : 


hg  -  (1.2)  f1/3 


(Eq.  2.9). 


2.  Compressive  Stress 

Several  empirical  relationshipa  have  also  been  developed  to 
estimate  peak  compreasive  stresses  developed  in  the  soil  as  a  function 
of  charge  weight  and  distance  from  the  detonation  point.  The  equation 
for  the  peak  compressive  stress.  opk  ,  in  MPa,  has  the  following  form: 

a.  -  K  (-8-)_n  (Eq.  2.10) 

Pk  W- 

where,  K  is  the  radial  distance  fro*  the  charge  in  meters,  W  is  the 
charge  weight  in  kilograms  per  delay,  K  is  a  ground  transmission 
constant  baaed  on  charge  confinement  conditions  and  local  geology,  and 
n  and  m  are  empirical  exponents.  When  m  is  taken  as  one  half,  the 
equation  is  referred  to  as  'square  root  scaling'  and  generally  matches 
results  from  row,  line  and  near  surface  charges.  When  m  is  taken  as 
one  third,  the  equation  is  referred  to  as  'cube-root  scaling'  and 
generally  matches  results  from  deeply  buried  point  charges  (Ambraseys 
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and  Hendron, 1968;  Dowding,1985;  Drake  and  Ingran,1981) .  In  the  field 
the  actnal  propagation  relationehipa  should  be  established  by  small 
scale  in  situ  explosive  tests. 


For  a  saturated  soil,  the  peak  compressive  stress,  o 
from  contained  point  charges,  is  given  by  (Lyakhov, 1961) : 


pk 


pk 


(58.9)  (- 


,1/3 


-1.05 


in  MPa, 


(Eq.  2.11). 


The  peak  compressive  stress,  ,  in  MPa,  from  concentrated 


charges  detonated  in  water,  is  given  by  (Cole, 1948): 


’pk 


(54.9)  (-flj) 


-1.13 


(Eq.  2.12). 


In  comparing  Equations  2.11  and  2.12,  it  can  be  seen  that  they 
differ  only  slightly  in  their  constants  which  indicates  that  the  peak 
compressive  stress  developed  in  a  saturated  soil  is  very  near  to  that 
developed  in  water.  However,  since  a  saturated  soil  is  a  two-phase 
medium,  different  behavior  can  be  expected  under  explosive  loadings. 

The  peak  compressive  stress,  ,  in  MPa,  from  concentrated 
charges  detonated  in  a  saturated  sandy  soils,  is  given  by  (Drake  and 
Little, 1983) : 

opk  -  (5.6  x  10-6)  (f)  (pt  Vc)  <^j73>'2’35  <Eq.  2.13) 

where,  f  is  a  ground  shock  coupling  factor  found  in  Figure  2.3  ,  p ^  is 
the  total  mass  density  in  kilograms  per  cubic  meter  and  V  is  the 
compressive  stress  wave  propagation  velocity  in  meters  per  second.  The 
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quantity  ' p  V  '  i*  commonly  referred  to  as  the  acoustic  impedance  or 
t  c 

specific  acoustic  resistance  of  the  medium. 

3.  Particle  Velocity 

A  widely  used  scaling  law  for  peak  particle  velocities,  V  ^  ,  in 

meters  per  second,  as  a  function  of  charge  weight  and  distance  from 
detonation  point,  has  the  following  form: 


V  ‘  1  ^ 


(Eq.  2.14) 


where,  R  is  the  radial  distance  from  the  charge  in  meters,  V  is  the 
charge  weight  in  kilograms  per  delay,  K  is  a  ground  transmission 
constant  based  on  charge  confinement  conditions  and  local  geology  and  n 
and  m  are  empirical  exponents.  When  m  is  taken  as  one  half,  the 
equation  is  referred  to  as  'square  root  sealing'  and  generally  matches 
results  from  row,  line  and  near  surface  charges.  When  m  is  taken  as 
one  third,  the  equation  is  referred  to  as  'cube-root  scaling'  and 
generally  matches  results  from  deeply  buried  point  charges  (Ambraseys 
and  Hendron.1968;  Dowding,1985;  Drake  and  Ingram, 1981) .  In  the  field 
the  actual  propagation  relationships  should  be  established  by  small 
scale  in  situ  explosive  tests. 

For  large  explosions  in  deep  soil  deposits,  Drake  and  Ingram  (1981) 
suggest  that  the  following  equations  be  used  to  estimate  the  peak 

per  second,  for  row  charge  explosive 

ar— 

configurations: 


particle  velocity,  V  ,  in  met< 

P* 


Pk 


(7.2)  (- 


,1/3 


-1.15 


(Eq.  2.15) 


where*  V  ie  the  eherge  weight  in  kilograas  per  delay  per  aster  of 
charge  row  length,  and  for  point  charge  explosive  sources: 


V  "  (8-0>  (p72)"2*3  (Eq-  2,16) 

where,  V  is  the  charge  weight  in  kilograas  per  delay. 

The  peak  particle  velocity,  Vpk  .  in  sieters  per  second,  froa  small 
buried  charges  in  saturated  sands,  is  given  by  (Long  et  si.,  1981): 

V  -  «>•«>  (Eq.  2-17>- 

Drake  and  Little  (1983)  recoaaend  the  following  equation  for  the 
peak  particle  velocity,  V*  in  asters  per  second,  for  contained  point 
charges  in  loose,  saturated  aande: 

vPk  "  (5'6)  <*>  (“iv3>~2*35  CBq.  2.18) 

where,  f  ie  a  ground  shock  coupling  factor  found  in  Figure  2.3. 

4.  mbmlifiXiffli 

Based  on  a  aeries  of  tests  using  buried,  contained,  point  charges 
in  loose,  saturated  sand,  Eok  (1977,1978b)  and  Studer  and  Kok  (1980) 
developed  an  eapirical  relationship  for  estiaating  residual  excess 
porewater  pressure  increases  in  the  field.  The  relationship  is 
expressed  in  nondiaensionalized  fora  as  a  function  of  charge  weight  and 
radial  distanoe  froa  the  charge: 
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1.65  +  0.64  In  ( 


.1/3 

V5 


(Eq.  2.19) 


where,  Lq  is  the  liquefaction  coefficient,  u^  is  the  residnsl  excess 

porewater  pressure,  a'  is  the  initial  effective  stress,  R  is  the  radial 

o 

distance  fron  the  charge  in  meters  and  I  is  the  charge  weight  in 
kilograms  per  delay.  The  relationship  given  by  Equation  2.19  is 
presented  in  Figure  2.4  for  various  porewater  pressure  ratios.  Field 
explosive  test  results  are  shown  in  Figure  2.5. 

G.  One-Dimensional  Stress  Wave  Props sat  ion 

Two  useful  equations  have  been  developed  from  the  theory  of 
elasticity  and  the  theory  of  stress  wave  propagation  in  a  linear 
elastic,  isotropic,  homogeneous  medium.  The  equations  were  derived 
considering  a  planar,  longitudinal,  compressive  stress  wave  propagating 
along  a  medium.  In  general,  these  equations  are  only  approximations 
for  soils  based  on  the  assumptions  in  their  derivations  and  can  only  be 
applied  for  small  strain  conditions.  However,  in  a  saturated  soil  with 
uadrained  loading  conditions,  the  response  is  primarily  governed  by  the 
water  phase,  making  the  assumptions  reasonable  and  justifying  the  use 
of  these  equations.  Rinehart  (1975)  indicates  that  these  relationships 
are  valid  on  the  wave  front  of  any  type  of  propagating  stress  wave,  but 
are  only  valid  for  planar  stress  waves  propagating  behind  a  wave  front. 


1.  Compressive  Stress 

A  linear  relationship  exists  between  the  instantaneous  stress  at 
any  point  in  an  elastic  medium  propagating  a  stress  wave  and  the 
particle  velocity  at  the  same  point.  The  equation  relating  the 
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CMpmiivt  stress.  «#  ,  sad  the  particle  velooity.  V  ,  is  given  by 
(Kolsky.1963;  Tiaoshenko  sad  Ooodier.1970) : 


<  P  Vc  >  Vp 


(Eq.  2.20) 


where,  p  is  the  aass  density  sad  is  the  conpressive  stress  wave 

propagstioa  velooity  through  the  aediaa.  The  qasatity  ' pV  *  is 

c 

coaaonly  referred  to  as  the  aconstic  iape dance  or  specific  acoustic 
resistance  of  the  aedina. 


The  coapressive  straia.  sq  .  developed  in  the  aedina  can  be 
deterained  fora  the  following  equation  (Kolsky.1963;  Tiaoshenko  and 


Goodier,1970) : 


(Eq.  2.21) 


It  can  be  seen  froa  Equation  2.21  that  the  coapressive  strain,  s  .  is 

c 

only  a  function  of  the  coapressive  stress  wsve  propagation  velocity 
through  the  awterial  and  the  particle  velocity  as  the  stress  wave 
propagates. 


h* 


% 
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III.  DEVELOPMENT  OF  A  NEW  EXPERIMENTAL  SOIL  DYNAMICS  FACILITY 


A.  Bilim  of  the  Experimental  Facility 

An  experimental  facility  was  designed  and  constructed  to  investigate 
the  transient  and  long-tern  porewater  pressure  response  of  saturated 
soils.  Sasiples  of  saturated  sand  can  be  prepared  at  varying  initial 
densities  and  effective  stresses.  The  experiaental  boundary  conditions 
are  for  one -dimensional,  confined*  conpressive  loadings  without 
drainage.  The  systea  is  capable  of  applying  conpressive  shock  loadings 
on  the  order  of  3S000  KPa  having  Billisecond  rise  tines  to  peak  stress. 

The  experiaental  facility  is  located  in  the  Departnent  of  Civil 
Engineering's  Geotechnical  Engineering  Research  Laboratory  at  Colorado 
State  University.  The  facility  is  divided  into  two  areas.  One 
contains  the  electronic  recording  equipaent  and  supporting  devices  and 
the  other  contains  the  experiaental  set-up  which  is  located  behind  a 
concrete  wall  for  safety.  Various  aspects  of  the  systea,  including  the 
design  considerations  and  construction,  will  be  discussed  in  the 
following  sections. 

B.  Cannon  for  ■Shock  Load  Application 

To  siaulate  the  conpressive  stress  conditions  associated  with 
explosive  loadings,  it  is  necessary  to  apply  a  high  aaplitude  stress 
having  a  Billisecond  rise  tine  to  peak  stress.  The  possibility  of 
using  a  coapressed  air  shock  loading  systea,  sinilar  to  that  used  in 
typical  shock  tube  experinents,  was  considered  but  then  abandoned  since 


it  did  not  Met  the  loading  requirements  of  this  experimental 
investigation.  Farther  investigations  into  the  literature  shoved  that 
direct  impact  methods  vonld  be  the  most  practical  way  to  obtain  the 
compressive  stress  pulse  desired,  vithont  using  explosives.  A  loading 
aystem  using  a  gas  powered  cannon  that  fired  a  6.80  cm  diameter,  10.15 
cm  long  stainless  steel  projectile  vas  designed  and  constructed.  The 
cannon  operates  using  compressed  nitrogen  gas  and  has  one  internal 
moving  part  called  a  'poppet  valve.'  The  principle  of  operation 
involves  applying  an  equal  pressure  on  both  sides  of  the  poppet  valve 
vhen  preparing  to  fire.  The  release  of  pressure  on  one  side  of  the 
valve  to  atmospheric  pressure  allows  a  rapid  release  of  the  remaining 
high  pressure  into  the  cannon  barrel  to  propel  the  projectile.  The 
projectile  velooity  was  varied  by  using  different  nitrogen  pressures  in 
the  cannon.  The  cannon  design  was  based  on  a  working  system  in  the 
laboratory  used  for  dynamic  loading  of  rock  core  specimens.  The  cannon 
vas  pressurized  and  fired  by  a  manually  controlled  mechanical  valving 
system.  Figure  3.1  shows  the  cannon. 

c.  CpnilnlM  Tube 

One  objective  of  these  experiments  was  to  examine  the  influence  of 
the  initial  effective  stress  on  the  sample's  porewater  pressure 
response.  A  controlled  method  of  applying  both  a  confining  pressure 
and  a  back  pressure  to  the  sample  was  developed.  Since  a  direct  impact 
method  of  stress  application  was  to  be  used,  consideration  was  also 


given  to  insuring  that  nearly  all  the  loading  was  applied  to  the  sample 
and  not  the  confining  tube  system. 


The  confining  tube  was  constructed  fro*  a  122  cm  long  section  of 
seamless*  stainless  steel  pipe  having  a  1.27  cm  wall  thickness  with  a 
8.90  cm  inside  diameter.  Welded  to  each  end  of  the  confining  tube  was 
a  30.48  cm  diameter.  2.34  cm  thick*  drilled*  stainless  steel  flange 
plate  used  to  bolt  the  entire  system  together.  High  strength  steel 
bolts  with  lock  washers  were  used  on  each  flange  plate.  The  confining 
tube  was  filled  with  distilled,  de-aired  water.  A  pressure  line  was 
attached  to  the  tube  at  the  midpoint  using  a  high  pressure*  one-way 
quick  disconnect  valve  in  series  with  a  high  pressure  plug  valve.  This 
type  of  connection  was  used  to  maintain  undrained  loading  conditions. 
Pipe  thread  fittings  were  need  in  connecting  the  valves  to  the 
confining  tube  for  a  high  pressure  seal.  The  confining  tube  was 
pressurised  by  a  small  hydraulic  pressure  vessel  filled  with  distilled, 
de-aired  water  and  connected  to  a  high  pressure  air  line.  The  pressure 
vessel  was  used  to  maintain  an  interface  between  the  air  pressure  and 
the  water.  The  distilled,  de-aired  water  in  the  confining  tube  and  the 
pressure  vessel  was  periodically  replaced.  Based  on  the  effective 
stress  required,  the  line  air  pressure  would  be  increased  to  the 
correct  confining  pressure  which  then  applied  pressure  to  the  confining 
tube  thro  .gh  the  interface  vessel.  The  interface  pressure  vessel  is 
shown  in  Figure  3.2. 

The  porewater  pressure  transducer  where  the  measured  pressure  would 
be  the  impact  stress  applied  to  the  sample.  The  transducer  was  located 
as  near  to'  the  sample  as  possible  in  the  confining  tube.  The 
transducer  was  flush  mounted  with  the  inside  of  the  confining  tube  wall 
and  a  high  pressure  seal  was  obtained  according  to  the  manufacturer's 


specifications 
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The  and  of  the  confining  tnbe  impacted  by  the  projectile,  was 
required  to  effectively  transmit  the  inpact  stress  to  the  confining 
pressure.  This  conld  not  be  accomplished  with  a  stiff,  rigidly  fixed 
loading  plate  sinoe  under  such  conditions  the  applied  stress  would  be 
transmitted  to  the  confining  tube.  In  meeting  this  requirement,  a 
stainless  steel  impact  piston  was  designed  that  was  not  directly 
attached  to  the  tube.  The  piston  was  allowed  to  move  under  the  applied 
load.  One  end  of  the  piston  was  machined  to  fit  within  the  end  of  the 
confining  tube  and  was  equipped  with  five  high  pressure  O-ring  seals. 
The  other  end  of  the  piston  was  machined  to  the  projectile  diameter 
with  a  bevelled  edge  for  placing  a  brass  loading  cap.  The  loading  cap 
was  used  to  prevent  excessive  deformations  of  the  exposed  end  of  the 
piston  from  repeated  impacts. 

A  2.54  cm  thick,  30.48  cm  diameter  stainless  steel  plate  was  bolted 
on  the  impact  end  to  hold  the  piston  in  place.  The  plate  was  fitted 
with  a  high  pressure  O-ring  seal  that  came  in  contact  with  the  end  face 
of  the  confining  tube.  The  hole  in  the  plate  for  the  piston  contained 
two  teflon  rings  to  minimize  friction  between  the  piston  and  the  end 
plate.  Details  of  the  impact  piston,  the  confining  tnbe  and  the 
porewater  pressure  transducer  location  are  shown  in  Figure  3.3. 

D-  Swlt  Cgptilatr 

An  important  requirement  in  the  design  of  the  sample  container  was 
I  to  have  undrained,  one-dimensional,  confined  compressive  loading 

conditions.  To  accomplish  this,  the  design  used  a  rigid  sample 
container  to  minimize  system  strain  effects  on  the  soil  sample.  It  was 

!  also  necessary  to  be  able  to  apply  a  back  pressure  within  the  sample 
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and  have  a  large  aspect  ratio  for  the  saaple  (i.e.  thin  saaple)  to 
ainiaize  any  inertia  effeots  that  aight  occur  with  the  passage  of  the 
stress  wave. 

The  saaple  container  was  constructed  froa  of  a  IS. 25  ca  long 
section  of  aeaaless,  stainless  steel  pipe  having  a  1.27  ca  wall 
thickness  with  an  8.90  ca  inside  disaeter.  Welded  to  each  end  was  a 
30.48  ca  diaaeter,  2.S4  ca  thick,  drilled,  stainless  steel  flange  plate 
used  to  bolt  the  entire  systea  together.  High  strength  steel  bolts 
with  look  washers  were  used  on  each  flange  plate.  Two  pressure  lines 
were  attached  to  the  container  at  the  aidpoint  using  a  high  pressure, 
one-way  quick  disconnect  valve  in  line  with  a  high  pressure  plug  valve 
on  each.  The  lines  were  looated  180  degrees  apart  and  used  in 
preparing  and  saturating  the  aaaple.  This  type  of  connection  was  used 
to  aaintain  undrained  loading  conditions.  Pipe  thread  fittings  were 
used  in  connecting  the  valves  to  the  saaple  container  for  a  high 
pressure  seal. 

In  pressurising  the  saaple  container,  it  was  necessary  to  prevent 
any  soil  grains  in  the  saaple  froa  aoving  out  into  the  two  pressure 
lines.  A  aodified  fitting  with  a  recessed  internal  slot  was  used.  A 
saall  porous  stone  was  sized  to  fit  snugly  in  the  recess  of  the  two 
fittings.  Bach  stone  was  securely  glued  in  place  with  a  silicone 
adhesive. 

Since  air  bubbles  in  the  saaple  oould  not  be  tolerated,  it  was 
necessary  to  flush  and  back  pressure  the  saaple  container  using 
distilled,  de-aired  water.  The  saaple  was  back  pressured  using  a  saall 
hydraulic  pressure  vessel  filled  with  distilled,  de-aired  water  and 
connected  to  a  high  pressure  air  line.  Figure  3.2  shows  the  pressure 
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vessel  which  was  tho  iu0  do  sign  as  that  used  for  tho  confining  tnbo. 
The  pro as arc  vessel  was  used  to  aaintain  an  interface  between  the  air 
pressure  and  the  water.  The  distilled,  de-aired  water  in  the  pressure 
vessel  was  periodically  replaced.  A  0.033  ca  thick  latex  rubber 
aeabrane  was  placed  between  the  saaple  container  and  the  confining  tube 
before  bolting  thea  together.  The  line  air  pressure  was  increased  to 
345  KPa  which  then  applied  the  back  pressure  to  the  saaple  container 
through  the  interface  vessel. 

The  porewater  pressure  transducer  was  located  where  an  average 
saaple  porewater  pressure  would  be  aeasured.  The  transducer  was  placed 
at  the  aidpoint  of  the  saaple  container.  The  transducer  was  flush 
aounted  with  the  inside  of  the  saaple  container  wall  and  a  high 
pressure  seal  was  obtained  according  to  the  manufacturer's 
specifications.  Although  the  transducer  was  not  located  in  the  center 
of  the  saaple,  results  from  its  position  would  be  indicative  of  the 
average  applied  coapressive  stress  amplitude.  The  porewater  pressure 
transducers  in  the  confining  tube  and  the  saaple  container  were  aligned 
and  positioned  on  the  saae  aide  and  in  the  saae  plane. 

The  end  of  the  saaple  container  that  attached  to  the  confining  tube 
was  aachined  saooth  to  seal  with  the  high  pressure  O-ring  on  the 
adjoining  face.  A  2.54  ca  thick,  30.48  ca  diaaeter  stainless  steel 
plate  was  attached  to  the  other  end  of  the  saaple  container  and  was 
fitted  with  two  concentric  high  pressure  O-ring  seals.  The  steel  plate 
was  bolted  securely  in  place  with  high  strength  steel  bolts  and  lock 
washers.  Details  of  the  saaple  container  and  porewater  pressure 
transducer  location  are  shown  in  Figure  3.4. 


57 


vil 


y) 


& 

*i 


To  effootivoly  apply  a  single  compressive  loading  to  the  soil, 
stress  wave  reflections  from  the  end  of  the  sasiple  container  needed  to 
be  miniaixed.  The  none n tun  trap  system  used  in  this  investigation 
typioally  reduced  reflected  stress  waves  to  less  than  about  50  percent 
of  the  applied  coapressive  stress. 

Based  on  these  observations  and  an  investigation  of  the  phenomenon 
of  stress  wave  propagation  through  aaterials,  several  criteria  were 
eatablished  to  aid  in  finding  an  efficient  aethod  of  stress  wave  energy 
disaipation: 


the  aaterial  should  have  an  acoustic  impedance  close  to  that  of 
water  or  the  soil-water  aixture.  The  better  the  acoustic  im¬ 
pedance  match  the  better  the  stress  wave  transmission  between 
mediuas. 


-  the  material  should  be  long  enough  to  capture  the  stress  wave. 


the  material  should  have  a  larger  contact  area  than  the  sphere, 
being  at  least  the  same  sixe  as  the  sample  cross  sectional  area 
to  be  the  most  effective,  and 


the  aaterial  should  be  free  to  aove  away  from  the  sample 
container  end  plate  as  the  stress  wave  passes  through  it.  This 
would  require  a  low  friction  support  mechanism  and  would  allow 
the  stress  wave  to  dissipate  its  energy  within  the  material 
apart  froa  the  saaple.  The  aaterial  also  must  not  be  allowed 
to  rebound  back  towards  its  original  position  where  it  could 
impart  a  stress  wave  due  to  its  own  motion. 


An  investigation  of  coaaercially  available  materials  to  satisfy  the 
first  requireaent  led  to  the  selection  of  a  solid  polyvinylchloride 
(PVC)  rod.  PVC  has  an  acoustic  impedance  just  above  that  of  water. 
Two  92  ca  long  sections  of  9  cm  diameter  PVC  rod  were  obtained  in 
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accordance  with  the  second  requirement.  To  satisfy  the  third 
requirement,  the  ends  of  the  rods  were  machined  and  trned  to  be  flush 
with  each  other  and  the  steel  end  plate. 

In  satisfying  the  fonrth  requirement,  a  low  friction,  free  moving 
support  system  was  designed.  A  305  cm  long  section  of  15  cm  diameter 
thin  walled  PVC  pipe  was  selected.  The  pipe  was  attaohed  to  a  305  cm 
long  section  of  8  cm  wide  steel  channel  seotion  which  was  fixed  to  and 
supported  by  two  heavy  duty,  adjustable  surveying  tripods.  The  solid 
PVC  rods  were  placed  inside  the  pipe  and  the  unit  was  aligned  with  the 
sample  container.  This  arrangement  provided  a  stationary  support 
system  that  allowed  the  rods  to  move  freely  as  the  compressive  stress 
wave  passed  through  them. 


To  simulate  the  type  of  loading  encountered  with  an  explosion. 


it  is  necessary  to  apply  a  compressive  stress  pulse  of  a  sufficient 
magnitude  with  milliseoond  rise  time  to  peak  stress.  The  measurement 
of  porewater  pressures  associated  with  this  type  of  loading  places 
speoific  requirements  on  the  transducers,  amplifiers  and  signal 
conditioners  necessary  to  accurately  monitor  the  pressure  changes  over 
the  time  domain  of  the  loading.  Considering  the  type  of  loading  to  be 
applied  in  the  laboratory  and  the  specific  requirements  of  this 
research  program,  a  number  of  criteria  were  established  to  evaluate 
coaauercislly  available  transducers  and  supporting  equipment.  The 
requirements  included  the  following: 
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the  transducer  nit  be  able  to  aeasure  the  porewater  pressure 
response  of  the  soil  saaple  independent  of  the  soil  skeleton 
stress* 


the  transducer  aust  have  a  linear  response  to  pressures  of  up 
to  35  MPa*  having  Billisecond  rise  times* 


the  transducer  aust  have  a  resonant  frequency  greater  than  the 
shock  facility  systea  where  it  is  used* 


the  transducer  aust  be  ale  to  survive  repeated  shock  loadings. 


the  transducer  aust  have  ainiaal  acceleration  sensitivity* 


r 

k 


the  transducer  aust  be  able  to  survive  soae  aaount  of  overpres¬ 
sure  and  still  have  a  linear  response* 


the  transducer  aust  able  to  survive  electronically  in  water  for 
at  least  several  days* 


R 


the  transducer  must  have  a  ainiaua  sensitivity  of  less  than 
seven  KPa  at  Full  Scale  Output  voltage, 

the  transducer  aust  be  able  to  measure  both  the  transient  and 
long-tera  porewater  pressure  responses* 

the  transducer  aust  be  able  to  consistently  reproduce  pressure 
responses  under  loading  conditions. 


the  transducer  aust  be  flush  aountable  in  the  steel  walls  of 
the  confining  tube  and  saaple  container,  having  a  high  pressure 
seal  of  at  least  17  MPa,  and 

the  transducer  aust  be  saall  in  sire  due  to  physical  space  lim¬ 
itations  and  soil  saaple  size. 


A  review  of  aany  instrumentation  companies’  products  led  to  the 
selection  of  a  porewater  pressure  transducer  manufactured  by  ENDEVCO, 
Inc.*  of  San  Juan  Capistrano,  California.  The  ENDEVCO  Model  8511a-5kNl 
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Mt  the  above  requirements.  It  it  a  piesoresist ive,  silicon  diaphragm, 
strain  gage  pressure  transducer.  The  signal  conditioner  units  were 
ENDEVCO  Model  4476.1a  and  the  power  supply  units  were  ENDEVCO  Model 
4470  for  these  transducers. 

The  transducer  was  modified  by  the  manufacturer  to  meet  the  first 
requirement.  A  one  millimeter  thick,  stainless  steel  circular  plate 
was  mounted  on  the  transducer  casing  oyer  the  diaphragm.  Hie  plate  had 
several  very  small  holes  (smaller  than  the  smallest  grain  size  of  the 
soil)  placed  in  it  in  a  star  pattern.  The  space  behind  the  plate  was 
filled  with  silicon  oil.  Sinoe  the  plate  was  not  in  contact  with  the 
transducer  diaphragm,  measured  stresses  would  correspond  to  the 
porewater  pressures  only.  Stresses  applied  by  the  soil  skeleton  would 
be  transmitted  through  the  points  of  contact  between  the  plate  and  the 
housing.  Figure  3.S  shows  the  modified  transducer.  Table  3.1  gives 
information  on  the  features  and  characteristics  of  the  transducer. 

2.  puiui  Igggjflat  Rtggrdgr 

To  simulate  the  type  of  loading  encountered  with  an  explosion, 
it  is  necessary  to  apply  a  compressive  stress  pulse  of  sufficient 
magnitude  with  millisecond  rise  time  to  peak  stress.  The  recording  of 
measured  data  associated  with  this  type  of  loading  places  specific 
requirements  on  the  instrumentation  necessary  to  accurately  monitor  the 
porewater  pressure  transducer  output  over  the  time  domain  of  the 
loading.  Considering  the  type  of  loading  to  be  applied  in  the 
laboratory  and  the  specific  requirements  of  this  research  effort,  a 
number  of  criteria  were  established  in  order  to  evaluate  cosusercially 
available  waveform  recorders.  The  requirements  included  the  following: 


Table  3.1  Manufacturer's  (ENDEVCO)  Specification  Data  Sheet  for 
the  Porewater  Pressure  Transducer 


The  Endevco  Model  151 1  it  a  rugged,  piezoresistive  pressure  transducer  for 
high  pressures.  It  has  a  H  in.  mounting  thread  and  ta  available  in  range*  from 
S  000  to50  000  psig  It  offera  high  shock  raaiatance.  high  overload  capability, 
and  low  photoaenattivity 

Endevco  pressure  tranaducera  feature  a  four-active  arm  strain  gage  bridge 
diffused  into  a  acuiptured  ailicon  diaphragm  for  maximum  sensitivity  and 
wideband  frequency  response  Self-contained  hybrid  temperature  compensa¬ 
tion  provides  stable  performance  over  the  wide  temperature  range  of  0  to  200*  F 
(-18  to  ♦83*CJ.  Endevco  tranaducera  also  feature  excellent  linearity  (even  to  3x 
range),  high  shock  resistance,  and  high  stability  during  temperature  transients. 


The  Model  851 1  is  widely  used  for  high  pressure  applications  such  as  studies  ! 
of  structural  loading  by  shock  waves  resulting  from  explosive  blasts,  pulsations  j 
in  hydraulic  systems,  and  combustion  problems  It  features  an  ablative  i 
coating  over  the  diaphragm  to  protect  against  particle  impingement 


SPECIFICATIONS  FOH  THE  8511  SEMES 

(Accoreme  to  ANSI  MV  ISA  SlanPWQa) 


PERFORMANCE 

•911-5* 

Ml  1-10* 

Ml  1-20* 

8911-MK 

RANGE  (pdg)  (t  pal  ■  «  M5  kP») 

5000 

10000 

30000 

50  000 

3ENSITIV1TV'  (mV/ptl  at  1000  Vdc) 

100  ±  030 

050  ±  015 

025  ±  008 

010  ±003 

RESONANCE  FREQUENCY*  (kHz) 

>500 

>500 

>500 

>500 

— 

0.6 

0.4 

05 

B 

08 

10 

Note* 

0.8 

02 

0.5 

05 

0.1 

05 

1.5 

12 

20 

100 

100 

100 

100 

ZERO  SHIFT 

0  1 

02 

Note4 

Note1 

0  1 

0.1 

0  1 

0  1  @  25  Ibf-  ft 

30 

3.0 

30 

30 

40 

40 

40 

40 

THERMAL  TRANSIENT  RESPONSE  (p»l/*F) 

(par  ISA-S37  10.  Para  0  7.  ProcaO  1) 

1.1 

1.5 

04 

0.1 

WARMUP  TIME  (aaconds)  1 

15 

15 

15 

15 

0  003 

0004 

0006 

0  003 

0001 

0002 

0  004 

0  003 

ksSII 

20  000 

30  000 

40000 

75  000 

'TtoufiWiacaMMMMat  1000  V0caacrtanon  It  may  bo  opafatatf  it  ottof  *011090  tKrt  mo  4n»n  voitapa  ihouKJ  So  tpacrfiad  ai  lima  ot  o*aar  Eryjwvco  Mo4«  **77  **?3  or  4470 
j  S*f*W  COAWMonor  1%  rocommarWto  M  Ito  oacrtSfton  aoutca 

<  ’MtfrovM  Mrtormonco  cm  to  owamoo  by  romovCI  Of  ttw  ftOUhvo  coating  o«i  ifto  OiapftroQm 
!  ■Oiornanpa  Mwtod  k>  7$  000  tot  A  atufi  m  cXiMilm  may  occu»  aiNH  ••poaura  10  praaauraa  atoova  M  000  pti 
*Owronea  ia  to  «0  000  pai  Zaro  afvn  at  7a  rang*  typteaxy  •  i«  ot  ?■  X  SO 

I  'TtfMw  spOCiAeaaoni  — noPto  0<t  aptoial  ortor 


NOT!  an  *a>uoa  a««  typical  o*a*  compynaatto  tampan  ttur#  rang*  netp  «*#*«  now* 


laca  at  mactwty  1 
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Table  3.1  (cont'd.)  Manufacturer's  (ENDEVCO)  Specification  Data 

Sheet  for  the  Porewater  Pressure  Transducer 


SPECIFICATIONS  FOR  MODEL  SS11 


»  *  *  (7S0  f  150)  LONG 
10  (2  4)  CMA  CABLE  s. 


_ L; 


MOO  *  025  THICK  J 
<14.2  OOH3I 
COPHR  PLATED 
STEEL  GASKET 


<1/0-24  UNP-2A  TMO  V~ 

COUNTERBORE 

37V  3*0  <9  52/9  M)  OIA 

>  030/060  (  7V1  52)  DEEP 


DIMENSIONS  IN 
INCHES  AND  (MILLIMETRES) 


ELECTRICAL 

ALL  MODELS 

EXCITATION 

Rated 

10  00  Vdc 

Maximum 

18  Vdc 

ELECTRICAL 

CONFIGURATIONS 

Four-active-arm  piezoresistive 
bridge 

POLARITY 

Positive  output  lor  increasing 
pressure 

RESISTANCE 

Input 

i  ooo  n 

Output 

isolation 

1  30011 

(ambient  conditions) 

>100  Mil  at  50  V 

PHYSICAL 

DEAO VOLUME 

0.0003  in.*  (0.0049  cmJ) 

MOUNTING 

Threaded  cate.  3/8-24  UNF-2A 
threaded  length  0  520  in. 

MOUNTING  TORQUE 

12  ±2  Ibfft  to  20  000  psi 

25  ±2  Ibf  *  ft  above  20  000  psi 

1  Ibfft  *  1  356N-m 

MATERIAL 

Caw 

Stainless  steel 

Exposed  to 

Measured  Fluid* 

Cate.  Paryiena  C.  epoxy, 
copper  gasket,  ablative 

CABLE 

Integral,  4-conductor,  shielded. 
Teflon-insulated  32  AWG 
(7/40)  silver  plated  copper 
conductors,  gray  PVC  jacket. 

30  in  (0.76  m)  long 

IDENTIFICATION 

851 1 -XX.  where  XX  denotes 
rated  pressure  in  psig 

WEIGHT 

0.4  oi.  (tig)  less  cable 

ENVIRONMENTAL 

TEMPERATURE  RANGE 

Compensated 

0*F  to  200° F  (-18*C  to  ♦93°C) 
Can  be  compensated  over  any 
200“ F  (1 10*C)  span  from  -65* F 
to  »225*F  on  special  order 

Maximum 

-65* F  to  ♦250°F  (-54“C  to 
♦121*C) 

VIBRATION 

1  000  g.  sinusoidal 

SHOCK 

HUMIDITY 
(per  MIL-STD-202E, 
Method  103B. 

20  000  g.  100  pS.  hall-sine 

Test  Condition  B) 

Isolation  resistanca  greater 
than  100  Mil  at  50  V.  External 
case  is  sealed  with  epoxy. 
Circuit  within  case,  vented 
through  tube,  is  coated  with 
Paryiena  C 

•Oea^ned  |<K  UN  mm  Dry  noo-coodiK.il  gaaea  Not  au>taw»  for  <iav  wtth  high  pH  or  low 

pH  hguKlt  long  term  •»potur«  to  water 
•poawa 

or  aipoaure  to  aolventt  which  may  attack 

Continued  product  Improvement  necessitates  that  Endevco  reserve  the  right  to  modify  these  specifications  without  notice 
RELIABILITY  Endevco  maintains  a  program  of  constant  surveillance  over  all  products  to  ensure  a  high  level  of  reliability  This 
program  includes  attention  to  reliability  factors  during  product  design,  the  support  of  stringent  Quality  Control  requirements,  and 
compulsory  corrective  action  procedures  These  measures,  together  with  conservative  specifications,  have  made  the  name  Endevco 
•ynonymouswith  reliability.  Endevco  s  Quality  and  Reliability  System  meets  the  requirementsof  MIL-0-9858A  and  MIL-STD-785A 
CALIBRATION.  Each  transducer  is  calibrated  at  room  temperature  tor  pressure  sensitivity  with  1000  Vdc  excitation  Special 
calibrations  available— consult  factory 

U  S.  Patent  Nos  4  065  970  and  4  093  933  apply  to  this  transducar 


the  wavefora  recorder  aust  be  able  to  interface  with  the 
computer  ayatea. 


the  waveform  recorder  anat  provide  digital  ontpnt  for  the 
pnrpoae  of  aaaa  storage  operationa  and  calculations. 


the  wavefora  recorder  anat  be  able  to  record  data  at  rates  op 
to  the  aegahertx  range,  with  variable  sampling  rates. 


the  wavefora  recorder  aust  be  capable  of  recording  enough  data 
pointa  to  represent  the  transducer  output. 


the  wavefora  recorder  aust  have  a  ainiaua  sensitivity  of  14  KPa 
when  in  line  with  the  transducer,  and 


-  the  wavefora  reoorder  aust  be  able  to  record  at  least  four 
channels  and  be  expandable  for  future  aodif icationa. 

A  coapariaon  of  aiailar  products  offered  by  various  coapaniea  led 
to  the  selection  of  a  digital  wavefora  reoorder  (GOULD/BIONATION  Model 
2803)  manufactured  by  GOULD,  Inc.,  Santa  Clara,  California  that  net 
these  requirements.  This  recorder  included  features  such  as  a  dual 
tiae  base  for  recording  at  two  different  sampling  rates,  the  ability  to 
select  the  nuaber  of  data  points  reoorded  by  each  time  base,  pretrigger 
and  delayed  node  recording,  and  the  ability  to  have  up  to  8  channels 
internally  trigger  and  record  at  the  same  time,  outputting  digital  data 
through  a  single  interface  port.  The  versatility  and  flexibility  of 
this  recording  system  proved  to  be  extremely  useful  during  the  conduct 
of  these  experiments.  Figure  3.6  shows  the  digital  waveform  reoorder. 
Table  3.2  gives  inforaation  on  the  features  and  characteristics  of  the 


wavefora  recorder 


Table  3.2  Manufacturer's  (GOULD /BIOMATION)  Specification 
Data  Sheet  for  the  Digital  Wave  Form  Recorder 


AMALOa  OUTPUTS:  MzptayVwtfaal(V): 

Diaplzy  Horizontal  (X): 

Driptay  Blanking  (Z,  2): 

Hot  Output: 

OmTALMnWACK:  Oota  Output 

nag  Output: 

u^WaaWa^n^i  Wa|H»»  ■ 

Output  Aaguoat: 

Output  Statu*: 

Baaoiri: 

Tima  Saaa  Control: 

SYSTMCQmBOIATWm  OianMh: 

Trigger: 

Output  Control: 

MSCSUAMOUS  Operating  Tamp.  Bongo: 

Potuar: 

Size: 

m  -  '  »  -■ 

wmgni. 

Warranty: 

Accessories: 


D/A  output  reproduces  complete  signal  in  memory  every 
2  ms.  0.8  V  fullscale.  adjustable  from  0.6  to  I  0  V  fullscale 
Vertical  position  control  allows  ±15  V  adjustment  of  base¬ 
line.  Vertical  XI  Expand  allows  full  scale  display  and  overlap 
of  multichannel  traces. 

Repetitive  2  ms  ramp  synchronized  with  V  output 
1  V  p-p  amplitude.  Horizontal  Expand  XI,  X2.  X5,  andXlO 
increases  ramp  slope  accordingly.  Maximum  amplitude 
clamped  at  ±  1 .5  V  Horizontal  position  control  allows 
±  1 .5  V  adjustment  of  origin 

Repetitive  4  fis  pulses  synchronized  with  X  ramp  retrace. 

Z  signal  is  nominally  a  0  lo  i  b  V  pulse  2  is  nominally  a 
♦  5  to  0  V  pulse  Also  usable  as  trigger  signals  to  synchronize 
the  internal  sweep  of  an  oscilloscope 

Front  panel  pushbutton  initiated  analog  output  whenever 
in  Dist^ay  mode  Duration  is  200  s  for  entire  memory  contents. 
0-1  V  p-p  amplitude,  adjustable  from  08  to  1 .2  V  p-p.  Reverts 
to  Display  output  after  single  plot  LED  "Ready"  indicator 
lit  when  unit  is  in  the  Plot  mode 


8  bits  parallel,  TTl.  levels  positive  true,  word  serial 
asynchronous  data  transfer  under  control  of  Flag  and 
Command  signals.  1  ms  to  500  MS/address.  2  ms  latency 
for  asynchronous  rates  slower  then  500  m s/point 
Positive  TTL  transition  indicates  data  word  on  output  lines 
can  be  read.  Minimum  pulse  width  0.2  ms. 

Negative  TTL  transition  requests  next  data  word.  Minimum 
pulse  width  0  1  ms 

Input  of  TTL  zero  (or  ground  connection}  requests  initiation 
of  digital  output  mode  LED  "Ready"  indicator  lit  when  unit 
is  in  Digital  output. 

Output  of  nonrwnaHy  f3  V  level  indicates  unit  is  in  the  digital 
output  mode 

Output  of  nominally  ±3  V  level  indicates  unit  is  recording. 

Ten  parallel  lines,  nine  lines  normally  TTL  high  when 
"EXT”  selected  on  front  panel  Time  Base  A  control  Input  of 
TTL  zero  (or  ground  connection)  on  any  one  line  selects  time 
base  in  conjunction  with  high  (open)  or  low  (ground)  on 
tenth  line  to  select  ms  or  ms  modifier 

Each  system  consists  of  a  single  master  control  unit  and  up 
to  three  slave  units  for  a  total  of  eight  channels  The  necessary 
rear  panel  cables  and  connectors  will  be  provided 
A  T  rigger  can  be  derived  from  any  of  the  eight  signals 
An  Output  Control  switch  (for  Digital  or  Plot)  is  provided 
concentric  with  each  vertical  position  control.  Pulling  the 
switch  out  disables  the  output  for  display.  Digital  or  Plot 
Only  one  channel  at  a  time  can  be  available  for  Digital  or 
Plot  outputs  AH  eight  channels  can  be  displayed  simul¬ 
taneously  for  comparison 


0-50*C 

100.  120. 220,  or  240  VAC  ±  10%.  approx  75  W.  50-400  Hz. 
Height.  5  25’ (13.4  cm)  Width.  12  75" (32  4  cm) 

Depth.  19"  (48  3  cm) 

Approx  25  lbs  (11  3  Kg) 

All  Biomation  products  are  warranted  against  defects  in 
materials  and  workmanship  for  one  year  from  date  of  delivery 
Each  unit  is  supplied  with  a  line  cord  and  a  copy  of  the 
Operating  and  Service  Manual 


Table  3.2  (cont'd.)  Manufacturer's  (GOULD/BIOMATION)  Speci¬ 
fication  Data  Sheet  for  the  Digital  Wave 
Form  Recorder 


1  MM//25  pF 

iuO  mV  io  50  V  selectabe  9  ranges  m  1-2-5  sequence 
100  V  peak.  35  V  rms 

Selectable  DC  or  AC  (low  frequency  cutoff  is  2  Hz) 
Selectable  Input.  Ground.  Record  or  Hold. 

Adjustable  ±  1  x  fullscale  input  sensitivity  range 

Two  LEDs  indicate  when  input  signal  level  exceeds  +  or  - 

the  full  scale  input  range 

DC  to  1  25  MHz  on  all  sample  rates. 


Beeofutfon:  8  bits  ( 1  part  m  256)  at  all  sample  rates 

Aperture:  The  Track-and-Hold  Circuit  exhibits  less  than  1  ns  aperture 

uncertainty. 

Maximum  Conversion  Rate:  5  MHz 

TIMHAMMOIOBY :  Sample  Inteoraf:  Internal  Single  or  dual  record  rates  with  independently 

selectable  sample  intervals  from  0.2  ns  to  100  ms  m  a  1  2-5 
sequence 

External  Selectable  for  input  of  sample  pulses  of  nominally 
♦  5  V  to  0  V  with  risetime  <  1 00  ns  Rate  continuously  variable 
from  0.2  ns  to  500 ms  per  address  11  delay  of  more  than 
500  fts  occurs  between  advance  pulses,  all  further  pulses 
will  have  2  ms  latency 

Memory  sba:  8  bits  x  2048  words  per  channel. 

Total  Record  Time:  2048  x  Sample  Interval,  vanes  from  04  ms  to  200  s  when 

using  Internal  sample  interval  selection. 

Recording  begins  at  (delayed)  Trigger  event  Recording 
stops  after  2048  samples  are  stored 
Recording  is  continuous  after  Arm  event  until  Trigger  is 
received,  which  starts  selected  delay  count.  End  of  delay 
stops  recording  process.  Delay  selection  allows  a  desired 
portion  of  the  data  recorded  prior  to  the  Trigger  event  to  be 
retained  in  memory 

Usable  in  either  Normal  or  Pretrigger  operation  with  Internal 
Time  Base  only  Time  Base  A  is  used  for  initial  part  of 
recording.  Switch  to  Time  Base  B  after  selected  Trigger  Delay 
Either  channel  may  be  held  while  the  other  channel  is 
updated. 

TfUQQOt  CMAVtACTftVUSTICS:  Trigger  Modes:  Auto  Recording  automatically  initiated  after  each  display 

output  sweep  (USEFUL  PRIMARILY  DURING  SETUP) 

Normal  Internal  or  External  trigger  event  accepted  du»>nq 
any  display  or  record  sweep  depending  on  Record  Mod  ? 
selected 

Single  1  rigger  event  accepted  only  after  Arm  function 

Trigger  Source:  Internal  T rigger  event  detected  from  either  signal  input 

External  Trigger  event  detected  from  I  font  panel  EXT  input, 
impedance  1  MI1//15  pF.  minimum  pulse  width  100  ns. 
minimum  amplitude  4U0  mV 

Trigger  Slope:  Selectable  ♦  or  - 

Trigger  Level:  Adiustable  *  1  x  fullscale  input  t  s  3  V  E  xternal) 

Trigger  Coupling :  Selectable  AC  or  DC 

Trigger  Deley:  Proportional  to  Sample  Interval  selected  or  input  Selectable 

via  three  front  panel  decade  switches  from  0  to  9990  sample 
intervals  in  increments  of  10  sample  intervals 

Arm  Function:  Initiated  bv  front  panel  momentary  pushbutton  or  via  rear 

panel  T7L  level  pulse  mpul  LED  "Ready'  indicator  lit  when 
unit  is  Armed  and  ready  10  be  triggered 


Duel  TlnwBwc: 

now. 


ANALOG  «WT  CHAAACTOttSTlCS 


Analog  Impodanca: 

Input  Voltaga  Ranga: 
Maximum  Input  Voltaga: 
Input  Coupling: 

Input  Connection: 

Input  Offset: 

Officsls  Indicators: 
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A  key  component  of  the  laboratory  instrumentation  system  is 
tke  computer.  Ike  computer  is  essential  for  controlling  the  electronic 
instrumentation  and  interpreting  the  recorded  data.  In  view  of  the 
importance  and  central  role  of  the  computer  in  integrating  the  overall 
laboratory  facility,  a  number  of  criteria  were  established  to  evaluate 
coamiercially  available  computer  systems.  These  criteria  included  the 
following: 


the  computer  system  must  be  able  to  interface  with  the  other 
laboratory  equipment,  serving  as  a  controller  where  necessary. 


the  computer  system  must  be  able  to  permanently  storing  the 
recorded  data  on  mass  storage  mediums. 


the  computer  system  must  be  able  to  analyze  and  graphically 
display  the  stored  data. 


the  computer  system  must  be  able  to  produce  hard  copy  output 
of  the  stored  data,  and 


the  computer  system  must  be  expandable  for  future 
modifications. 


Based  on  a  comparison  of  the  products  offered  by  several  companies, 
HEWLETT-PACKARD  was  chosen  to  supply  the  computer  system  as  it  best  met 
the  above  requirements.  The  computer  system  originally  obtained  at  the 
start  of  this  research  effort  consisted  of  a  desktop  computer.  Model 
HP-983 5A  with  CRT  display  and  an  eight  pen  plotter.  Model  HP-9872C. 
Ike  system  interface  directly  with  the  digital  waveform  recorder  and 
the  university's  main  computer  facility.  Since  the  start  of  this 
research  effort,  a  dot-matrix  line  printer.  Model  HP-82905B  and  a 
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floppy  disc  Mt*  storage  unit.  Model  HP-9885M  have  been  added.  This 
system  has  proved  to  be  flexible  and  capable  of  handling  the  needs  of 
the  laboratory  facility.  Figure  3.6  shows  the  computer  and  supporting 
peripherals  devices. 

4.  Triaaerina  Mechanism 

An  external  triggering  aechanisa  was  used  to  initiate  the 
recording  of  transducer  output  data  by  the  digital  waveform  recorder. 
An  HP-5300A/5304A  digital  tiner/counter  unit  was  used  to  control 
triggering  and  determine  the  projectile  impact  velocity.  Just  before 
the  projectile  iaqiacted  the  loading  cap.  it  opened  an  electrical 
circuit  that  started  the  tiaer/counter.  Vhen  the  projectile  impacted 
the  brass  loading  cap,  the  circuit  was  closed,  the  timer/counter 
stopped  and  a  dc  triggering  voltage  was  sent  to  the  digital  waveform 
recorder.  The  triggering  system  initiated  data  acquisition  immediately 
upon  the  application  of  the  compressive  loading  to  the  piston. 


IV.  EXPERIMENTAL  INVESTIGATION 


A.  Description  of  Monterey  No.  0/30  Sand 

1.  jam  of  JUlixlai 

The  soil  need  in  this  investigation  was  obtained  from  Lone  Star 
Indnstries  Inc . ,  Oakland.  California  in  1982  by  the  United  States 
Department  of  the  Interior.  Bnrean  of  Reclamation  in  Denver.  Colorado 
and  Colorado  State  University.  A  quantity  of  900  kilograms  was 
aequired  with  22S  kilograms  being  stored  at  Colorado  State  University 
for  these  experiments.  The  soil  is  a  clean,  uniform,  subangular  fine 
sand  processed  from  material  dredged  from  the  beach  north  of  Monterey, 
California  and  is  used  in  sand  blasting  operations.  The  company 
designation  for  the  sand  is  'Monterey  No.  0/30'  which  is  based  on 
gradation  requirements  and  this  term  will  be  used  throughout  the  text 
in  reference  to  the  sand.  A  detailed  description  of  the  sand  and  its 
behavior  under  cyclic  loading  is  given  by  Muxxy  (1983)  and  Charlie  et 
al.  (1984). 

2.  Pi&Rtxntlga  at  tin  Bvlfc  fluml* 

Five  43  kilogram  bags  were  chosen  at  random  from  the  lot  and 
transported  from  the  Bureau  of  Reclamation  to  the  Department  of  Civil 
Engineering's  Geotechnical  Engineering  Research  Laboratory  at  Colorado 
State  University.  The  sand  was  thoroughly  mixed  to  insure  a  uniform 
distribution  of  material  so  that  samples  used  in  the  experiments  would 
be  representative  of  the  bulk  quantity.  After  mixing,  the  sand  was 
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divided  into  23  kilogram  quantities  and  stored  in  oanvas  bags  with 
plastic  liners.  Randoa  samples  were  chosen  to  be  used  in  analyzing  the 
physical  properties  of  the  sand. 

3.  Physical  Properties  of  Monterey  No.  0/30  Sand 

Several  tests  were  performed  to  investigate  various  physical 
properties  of  the  Monterey  No.  0/30  sand.  All  tests  were  conducted 
according  to  the  American  Society  for  Testing  and  Materials  (ASTM) 
standard  laboratory  procedures  where  applicable  and  include  the 
following: 

-  grain  size  analysis  (ASTM  D422) . 

-  classification  (Unified  Soil  Classification  System.  ASTM  D2487). 

-  specific  gravity  (ASTM  D8S4) , 

-  relative  density  determination  (ASTM  D2049) . 

-  photomicrograph  (Bureau  of  Reclamation), 

-  spectrographic  analysis  (Bureau  of  Reclamation)  and 

-  skeleton  stress-strain  curves  (Hendron,1963;  Whitman  et  al.  1964). 

Selected  test  results  and  analysis  are  given  in  Chapter  V. 

B.  Variation  of  Parameters 

In  considering  the  number  of  factors  observed  to  influence 
liquefaction  presented  in  Chapter  II,  it  was  necessary  to  limit 
variables  to  be  investigated  during  this  research  effort.  Three 
important  and  significant  parameters  that  could  be  varied  in  a 
controlled  manner  in  the  laboratory  were  selected  for  study: 
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-  the  initial  pecked  relative  density  of  air  dry  samples. 

-  the  initial  effective  stress  on  the  samples.  and 

-  the  number  and  intensity  of  applied  compressive  shock  loads. 

The  initial  relative  densities  of  the  samples  used  in  this 
experimental  investigation  were  approximately  0,  20,  40,  60  and  80 

percent.  Each  of  these  five  relative  densities  vere  subjected  to  four 
different  initial  effective  stresses,  which  were  86  KPa,  172  KPa,  345 
KPa  and  690  KPa.  For  eaoh  sample  the  back  pressure  was  maintained  at 
345  KPa  with  the  confining  pressure  adjusted  according  to  the  required 
effective  stress.  The  first  series  of  samples  was  loaded  by  four 
'high'  stress  impacts  (between  2.6  MPa  and  8.5  MPa  per  impact)  and  the 
second  series  of  samples  was  loaded  by  six  'low'  stress  impacts 
(between  0.10  MPa  and  4.2  MPa  per  impact). 

c.  fliaali  Pmmtioii 

The  preparation  of  each  sample  used  in  the  experiments  followed  a 
prescribed  method  of  placement  to  insure  uniformity  and  consistency 
throughout  the  course  of  the  research.  Each  sample  was  rsndomly  taken 
from  the  23  kilogram  bags  in  which  the  bulk  of  the  material  was  stored. 
All  samples  were  placed  in  the  sample  container  at  the  air  dry  moisture 
content  of  less  than  0.10  percent. 

i*  ElMuavkl  fothgdt 

In  preparing  samples,  two  methods  of  placing  the  sand  were  used 
to  meet  the  relative  density  requirements.  One  method  used  a  funnel 
while  the  other  involved  a  combination  of  funnel  placement  and 


controlled*  layered  compaction.  Before  placing  sand  in  the  sample 
container,  carbon  dioxide  gas  was  flushed  through  both  of  the  pressure 
valving  systems  and  the  sample  container.  Then  the  placement  of  the 
sand  began  with  the  carbon  dioxide  gas  continually  flowing  until  the 
sample  was  completed.  The  vacuum  line  was  attached  to  the  remaining 
inlet  valve  which  was  kept  closed  until  the  sample  placement  was 
completed  and  the  membrane  was  in  place. 

a.  Emafll  Pi  **«■««* 

To  obtain  a  very  low  density  sample,  near  zero  relative 
density,  the  sand  was  placed  using  an  18  cm  diameter  metal  funnel  with 
a  1.27  cm  diameter  by  23  cm  long  spout.  The  larger  part  of  the  funnel 
held  about  one  half  the  amount  of  sand  required  to  fill  the  sample 
container. 

The  empty  funnel  was  placed  in  the  sample  container  with  the  bottom 
of  the  spout  held  against  the  container  bottom.  Then  about  half  of  the 
sand  was  poured  in.  The  funnel  was  slowly  lifted  upward  and 
simultaneously  rotated  inside  the  sample  container  while  keeping  the 
spout  about  1.27  cm  or  less  above  the  plaoed  material.  During  this 
motion  additional  sand  was  placed  in  the  funnel  for  continuous  sample 
placement. 

b.  SadBiflunifiilfla 

The  preparation  of  samples  at  relative  densities  greater  than 
zero  percent  was  accomplished  using  a  method  presented  by  Ladd  (1978) 
called  weight.  Each  layer  was  placed  using  the  funnel  approach 
described  previously.  The  layers  were  individually  compacted  to 


successively  higher  percentages  of  the  final  sample  density,  varying 
linearly  by  layer.  Vhen  the  final  layer  was  placed*  the  entire  sample 
was  at  the  required  density.  Based  on  Ladd's  (1978)  findings  it,  was 
decided  that  ten  layers  of  soil  would  be  used  and  that  the  first  layer 
would  be  five  percent  underconpacted.  The  unde rcompact ion  nethod 
provided  a  uniform  density  across  the  extent  of  the  saatple.  The  tools 
and  equipment  used  are  shown  in  Figure  4.1. 

2.  ftmfellM  Rli&mai 

Vhen  all  the  sand  had  been  poured*  the  top  of  the  saaiple  was 
leveled.  Loose  sand  grains  on  the  flange  plate  were  placed  into  the 
sample  container  and  the  top  releveled.  Then  a  piece  of  0.00254  cm 
thick  aluminum  foil  about  10  cm  in  diameter  was  centered  over  the 
sample  and  held  in  place  with  a  small  amount  of  vacuum  grease  where  it 
contacted  the  flange  plate.  The  aluminum  foil  was  used  to  provide  a 
buffer  between  the  sand  and  the  membrane  reducing  membrane  penetration 
effects  and  protecting  the  membrane  from  puncture.  After  the  foil  was 
in  position*  a  thin  layer  of  vaouum  grease  was  placed  on  the  flange 
plate  where  the  rubber  membrane  would  be  located.  The  vacuum  grease 
was  evenly  distributed  and  then  a  16  cm  diameter,  0.033  cm  thick 
circular  latex  rubber  membrane  was  centered  over  the  sample.  The 
membrane  was  flattened  by  carefully  working  around  its  circumference 
until  there  was  about  1.27  cm  left  to  seal.  A  final  check  for  sand 
particles  on  the  flange  plate  was  made  and  the  last  section  of  the 
membrane  was  pressed  down  and  sealed.  The  vacuum  line  was  then  opened 
to  draw  down  the  membrane  into  final  position  while  the  carbon  dioxide 
gas  continued  to  flow  through  the  ssmple  container.  With  the  .temple 


under  vicuna  (35  KPa),  several  aeasurenents  across  the  aeabrane  were 


aade  to  deteraine  the  placed  dry  density  of  the  sample. 

D.  System  Assembly 

Once  the  saaple  container  was  prepared,  the  next  step  was  to- 
asseable  the  entire  system  as  a  unit.  The  O-ring  seal  on  the  confining 
tube  flange  plate  was  lightly  coated  with  vacuua  grease  and  the  tube 
was  filled  with  distilled,  de-aired  water.  The  sample  container  was 
then  lifted  using  a  chain  falls  suspended  from  the  overhead  roller 
support  systea  to  a  point  above  the  confining  tube,  then  aligned  and 
slowly  lowered.  When  the  saaple  container  was  about  2.54  cm  above  the 
flange  plate,  the  space  between  the  flange  and  aeabrane  was  filled  with 
distilled,  de-aired  water  and  then  the  saaple  container  was  lowered  to 
its  final  position.  Before  bolting  and  sealing  the  flange  plate 
connection,  an  alignment  check  of  flange  plate  bolt  holes  was  aade. 
The  bolts  were  tightened  in  a  specified  pattern  until  the  lock  washers 
were  moderately  coapressed.  Then  the  bolts  were  tightened  to  their 
final  position.  These  steps  were  taken  so  that  the  high  pressure  0- 
ring  seal  would  seat  evenly  about  the  flange  plate. 

With  the  saaple  container  and  confining  tube  bolted  together, 
another  chain  falls  was  attached  to  the  opposite  end  of  the  confining 
tube.  By  alternately  raising  and  lowering  the  two  ends  of  the  unit, 
the  entire  systea  aligned  horizontally  with  the  cannon  barrel.  Two 
pipe  hangers,  suspended  from  the  overhead  roller  support  system,  were 
positioned  at  each  end  of  the  assembly.  The  hangers  were  attached  to 
the  confining  tube  and  both  chain  falls  removed.  The  confining  tube 
and  saaple  container  were  then  supported  by  the  pipe  hangers  and  free 


■oving  roller  support  system.  The  assembled  system,  ready  for  loading, 
is  shown  in  Fignres  4.2  and  4.3. 


A  pressure  line  was  attached  to  the  confining  tube,  35  KPa  pressure 
was  introduced  into  the  line,  and  the  plug  valve  opened.  The  carbon 
dioxide  pressure  line  was  connected  to  the  bottom  of  the  sample 
container  and  a  line  open  to  the  atmosphere  was  connected  to  the  top  of 
the  sample  container.  The  valves  on  both  sides  were  closed  at  this 
time.  The  carbon  dioxide  line  was  opened  under  a  controlled  flow  rate 
at  a  low  pressure  for  about  15  seconds.  Then  the  line  connected  to  the 
top  of  the  sample  container  was  opened  to  the  atmosphere  and  the  carbon 
dioxide  gas  continued  to  flow  across  the  sample  for  30  minutes. 

After  30  minutes,  the  carbon  dioxide  pressure  line  was  removed  and 
the  plug  valve  closed.  The  valving  system  attached  to  the  interface 
pressure  vessel  was  then  purged  of  any  air  bubbles  using  a  distilled, 
de-aired  water  supply  suspended  about  70  cm  above  the  sample  container. 
Then  the  bach  preasure  line  was  connected  to  the  bottom  of  the  sample 
container,  the  plug  valve  opened  and  distilled,  de-aired  water  was 
flushed  through  the  sample.  The  line  coming  out  of  the  top  of  the 
sample  container  was  put  into  a  1000  milliliter  graduated  cylinder  to 
collect  the  outflow.  After  about  three  pore  volumes  of  distilled,  de- 
aired  water  had  passed  through  the  sample,  the  top  and  bottom  plug 
valves  were  closed  and  the  top  line  removed.  The  confining  pressure 
and  the  back  pressure  were  incrementally  increased  to  obtain  the 
required  effective  stress.  The  sample  was  allowed  to  saturate  and 
consolidate  for  a  minimum  of  6  hours. 
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Figure  4.2  Experimental  Shock  Facility 


Prepared  for  Loading  (View 

from  r.nnnnnl 


Once  the  saturation  period  was  coapletad>  it  was  then  necessary  to 
cheek  the  sample's  porewater  pressure  response  and  determine  the  degree 
of  saturation.  The  porewater  pressure  transducers  were  connected  to 
the  signal  conditioner  panel  and  the  excitation  voltage,  amplifier  zero 
and  bridge  balance  were  checked  and  adjusted.  The  back  pressure  was 
increased  about  35  KPa  and  the  inlet  valve  closed.  An  output  voltage 
was  reoorded  for  each  transducer. 

The  confining  pressure  was  incrementally  increased  to  the  required 
effective  stress,  and  the  output  voltage  of  each  transducer  was  noted. 
The  degree  of  saturation  was  determined  by  comparing  the  sample's 
porewater  pressure  response  to  the  increase  in  oonfining  pressure.  If 
the  ratio  of  these  two  values  was  unity,  corresponding  to  a  100  percent 
porewater  pressure  response,  then  the  sample  was  considered  saturated. 
In  this  investigation  it  was  not  always  possible  to  obtain  a  porewater 
pressure  ratio  equal  to  one  even  though  the  same  saturation  procedure 
was  consistently  used  for  each  sample.  However,  measured  compressive 
stress  wave  propagation  velocities  through  the  samples  showed  they  were 
saturated.  It  was  observed  that  the  porewater  pressure  response  varied 
in  a  predictable  manner  with  increasing  sample  density  and  effective 
stress.  A  presentation  and  discussion  of  this  observation  is  given  in 
Chapter  V. 
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When  the  confining  pressure  and  the  back  pressure  were  restored  to 
their  desired  values,  the  two  plug  valves  (one  on  the  confining  tube 
and  one  on  the  sample  container)  were  closed  to  prevent  drainage 
through  the  pressure  lines.  The  valves  were  left  closed  throughout  the 
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loading  itqwnct  for  each  experimental  ran,  to  maintain  undrained 
loading  conditions. 

The  confining  tube  and  saaple  container  assenbly  were  aligned  with 
the  cannon  barrel  and  leveled.  The  brass  loading  cap  was  placed  on  the 
iapact  piston.  The  triggering  aechanisa  was  set.  its  operation  checked 
and  the  tiaer/eonnter  was  reset.  The  PVC  rods  in  the  aoaentna  trap 
asseably  were  aligned  with  and  oentered  on  the  steel  end  plate  on  the 
saaple  container  and  the  tripods  leveled.  A  final  check  of  the  entire 
systea  alignaent  was  aade. 

The  output  of  the  two  transducers  was  adjusted  to  a  value  of  zero 
volts  to  gain  aaxiaua  benefit  of  the  digital  waveform  recorder's  range. 
The  controls  on  the  wavefora  recorder  were  set  and  adjusted  and  then 
verified  by  the  computer.  The  computer  then  initiated  a  prograa 
sequence  to  handle  the  data  transfer  froa  the  digital  wavefora  recorder 
for  each  iapact  and  store  the  output  on  a  mass  storage  aediua.  A  final 
check  of  all  systea  coaponents  was  made  before  the  loading  phase. 

The  next  step  was  to  load  the  saaple.  Vith  the  projectile  inserted 
in  the  cannon  barrel,  the  cannon  was  pressurised  and  the  reservoir 
pressure  was  verified  with  a  gage  separate  froa  the  nitrogen  supply 
bottle.  Then  the  pressure  behind  the  poppet  valve  was  reduced  to  the 
ataospheric,  which  opened  the  poppet  valve,  releasing  the  high 
pressure  into  the  cannon  barrel  to  propel  the  projectile.  The  value  on 
the  tiaer/counter  was  recorded  and  the  computer  initiated  the  data 
transfer  and  storage.  The  tine  recorded  froa  the  tiaer/counter  was 
used  to  deteraine  the  iapact  velocity  of  the  projectile.  This  sequence 
was  repeated  for  each  iapact  loading. 


Vhia  u  experiaental  ran  «u  conpleted,  the  entire  systen  wee 
disasseabled  in  reverse  order.  A  viobiw  wee  applied  to  the  staple 
container  to  drain  the  pore  fluid.  Vith  the  staple  container  suspended 
froa  the  chain  falls,  the  aeabrane  and  foil  paper  were  removed.  the 
soil  staple  reaoved,  and  the  staple  container  thoroughly  cleaned  and 
dried.  Each  of  the  porous  stones  were  subjected  to  a  vacuua  to  reaove 
water.  The  next  staple  was  prepared  and  loaded  according  to  the 
procedure  previously  described. 


V.  EXPERIMENTAL  RESULTS 


A.  Physiol  Properties  of  Monf  rev  No.  0/30  Sand 

Basic  physical  sad  index  properties  for  Monterey  No.  0/30  sand  were 
evaluated  according  to  accepted  standard  laboratory  testing  procedures 
set  forth  by  the  Aaierican  Sooiety  for  Testing  and  Materials  (AS1M) .  A 
summary  of  the  tests  perforated  is  presented  given  in  Table  5.1. 
Details  of  the  grain  sixe  analysis*  classification  and  relative  density 
tests  are  given  by  Mnxzy  (1983)  and  Charlie  et  al.  (1984). 

1.  Grain  Sixe  Analysis 

A  grain  size  distribution  analysis  (ASTM  D422)  was  perforated. 
The  results  obtained  indicate  a  fine,  uniform,  poorly  graded  sand  with 
less  than  one  percent  of  the  material  being  finer  than  0.150  mm  in 
size.  The  D^  particle  size  for  this  material  is  0.45  mm.  A  grain 
size  curve  showing  the  distribution  of  the  particle  sizes  is  shown  in 
Figure  5.1. 

2.  tlUAiflCtUPB 

From  the  results  of  the  grain  size  analysis  the  uniformity 

coefficient*  Cq  ,  was  determined  to  be  1.65  snd  the  coefficient  of 

curvature*  C  ,  was  determined  to  be  1.00.  Based  on  these  results  and 
c 

the  grain  size  curve*  the  soil  was  classified  as  sn  SP  material  (sand* 
poorly  graded  with  little  or  no  fines)  according  to  the  Unified  Soil 
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Table  5.1  Physical  Properties  of  Monterey  No.  0/30  Send 
(Muzzy, 1983;  Charlie  et  si., 1984) 


DSCS  Classification 

SP 

Specific  Gravity 

2.65 

Particle  Size  Data: 

D10 

0.29  mm 

°30 

0.38  mm 

D50 

0.45  mm 

cu  (1) 

1.65 

Cc  (2) 

1.00 

%  Passing  #100  Sieve  (3) 

0.05  % 

Relative  Density  Test  Data: 

Dry  Dnit  Weight: 

Maximum 

1700  Kg/M3 

Minimum 

1470  Kg/M3 

Void  Ratio: 

Maximum 

0.803 

Minimum 

0.563 
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Classification  System  (ASTM  2487).  From  the  results  of  a  specific 
gravity  test  (ASTM  D8S4) ,  a  value  of  2.65  was  determined  for  the  sand. 


A  relative  density  test  (ASTM  D2049)  was  performed  to  determine 
the  minimum  and  maximum  dry  density  limits  for  the  sand.  This  test  was 
performed  at  the  Bureau  of  Reclamation  Soil  Mechanics  Laboratory  in 
Denver,  Colorado,  using  a  calibrated  shaker  table.  As  is  shown  in 
Table  5.1,  the  range  of  dry  densities  varied  by  230  Kg/M^. 
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Figure  5.2  shows  a  photomicrograph  of  the  sand  made  with  a 
scanning  electron  microscope  at  the  Bnrean  of  Reclamation  Soil 
Mechanics  Laboratory  in  Denver,  Colorado.  The  soil  grains  are  uniform 
in  sixe  and  are  subangular  to  subrounded  in  shape.  Figure  5.3  shows  a 
spectrographio  analysia  performed  at  the  same  laboratory  facility  to 
investigate  the  mineralogic  content  of  the  sand.  The  predominant 
mineral  constituent  present  is  silicon  with  several  others  existing  in 
substantially  smaller  amounts. 


Static  one-dimensional,  confined,  compression  tests 
(Hendron,1963;  Whitman  et  al.,1964)  were  performed  on  sir  dry  samples 
of  Monterey  No.  0/30  sand.  The  results  of  tests  were  used  to  obtain 
stress-strain  information  for  the  soil  skeleton  and  to  determine  the 
constrained  modulus  of  the  skeleton  to  be  used  in  the  finite  difference 


approximation  presented  in  Chapter  VI.  To  simulate  the  initial  stress 
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Figure  5.3  Spectrographic  Analysis  Results 

for  Monterey  No.  0/30  Sand  (Muzzy, 
1983;  Charlie  et  al.,1984) 


conditions  in  the  experimental  investigation,  a  corresponding  initial 
effective  stress  was  applied  to  each  sample.  The  compressive  strain 
values  were  referenced  to  the  initial  applied  stress.  Each  sample  was 


loaded  and  unloaded  in  increments  two  times  to  develop  the  stress- 
strain  relationship  for  the  soil  skeleton.  From  these  results,  a 
constrained  modulus  for  loading  and  unloading  was  determined. 

Two  sets  of  tests  were  run  for  comparative  purposes.  Samples 
tested  had  relative  densities  of  40  and  80  percent.  Each  sample  waa 
tested  at  initial  effective  stresses  of  86  KPa  and  690  KPa.  The  data 
obtained  for  the  relative  densities  and  effective  stresses  examined 
provides  practical  limits  for  analysis.  The  skeleton  stress-strain 
curves  for  Monterey  No.  0/30  sand  are  shown  in  Figures  5.4,  5.5,  5.6 
and  5.7  and  a  summary  of  the  numerical  results  is  given  in  Table  5.2. 
The  results  show  that  the  soil  skeleton  stiffness  increases  with 
increasing  initial  effective  stress  and  density.  Hysteresis  between 
the  loading  and  unloading  curves  decreases  with  increasing  initial 
effective  stress  and  density. 

b.  S.tiUc  Rgipgaiy 

Before  loading  of  each  sample  the  porewater  pressure  response  was 
checked  to  determine  the  degree  of  saturation.  This  was  done  by 
increasing  the  confining  pressure  on  the  sample  and  monitoring  the 
sample's  porewater  pressure  response  without  drainage.  The  ratio  of 
the  sample  porewater  pressure  response  to  the  increase  in  confining 
pressure  is  termed  the  ' C-parameter'  (Lambe  and  Whitman,  1969)  for  a 
one-dimensional  confined,  compressive  loading  of  a  saturated  soil  with 
undrained  conditions.  A  ratio  of  one  indicates  a  saturated  sample  and 
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Table  5.2  Constrained  Loading  and  Unloading  Modulus  for  Static. 

One-Diaensional,  Confined  Compression  Tests  on  Dry 
Monterey  No.  0/30  Sand 


D 

a 

D. (1> 

D  (2) 

r 

o 

1 

n 

(%> 

(IPs) 

(KPa) 

(KPa) 

40 

172 

61500 

307500 

40 

690 

156700 

617000 

80 

172 

86190 

431000 

80 

690 

191700 

958500 

Note:  (1)  Constrained  nodnlos  for  loading 

(2)  Constrained  siodnlna  for  unloading  taken  as  5  times 

constrained  modulus  for  loading  in  finite  difference 
analysis 


\ 

1 

I 


value*  less  than  one  Indicates  that  the  saaple  is  not  saturated  or  has 
a  very  stiff  soil  skeleton.  For  an  initial  effective  stress  of  86  KPa. 
a  ratio  of  one  vas  consistently  obtained.  However,  a  C-parameter  of 
one  vas  not  obtained  for  higher  initial  effective  stresses.  Since  the 
preparation  and  saturation  process  was  identical  for  each  saaple,  it 
vas  assuned  that  the  porevater  pressure  ratios  obtained  indicated  of  a 
saturated  condition.  An  exaaination  of  the  compressive  stress  wave 
propagation  velocity  through  saaples  verified  this  assuaption.  The 
calculated  coapressive  stress  wave  velocities  were  close  to  1500  aeters 
per  second  in  all  saaples  investigated.  This  is  the  value  that  would 
be  expected  for  saturated  conditions. 

Throughout  the  experiaental  investigation  it  was  noted  that  the 
porevater  pressure  ratio  varied  in  a  predictable  aanner  vith  variations 
in  effective  stress  and  relative  density  (Figure  5.8).  In  considering 
this  observation  and  those  previously  discussed,  it  is  believed  that 
the  porevater  pressure  ratio  response  noted  can  be  attributed  to 
changes  in  the  soil  skeleton  stiffness  vhich  increases  vith  increasing 
initial  relative  density  and  effective  stress.  Accordingly,  all 
saaples  vere  considered  to  be  saturated. 

C. 

The  pressure-tine  histories  represent  the  porevater  pressure 
transducer  responses  to  applied  shock  loadings  as  function  of  tiae. 
They  include  both  the  peak  and  long-term  response  for  the  confining 
pressure  and  the  saaple  porevater  pressure.  Selected  pressure-time 
histories,  representative  of  the  behavior  observed  in  this  experiaental 
investigation,  are  shovn  in  Figures  5.9  through  5.28.  A  sunary  of 


nuaerical  results  f  roa  the  pressure-tine  histories  of  ill  saaple s 
investigated  is  given  in  Chapter  VI,  Tables  6.2  through  6.12. 

Figures  5.9  through  5.18  show  the  experiaental  results  for  the 
'40%'  relative  density  series  at  an  effective  stress  of  172  KPa. 
Figures  5.19  through  5.28  show  the  experiaental  results  for  the  '80%' 
relative  density  series  at  an  effective  stress  of  172  KPa.  The 
confining  pressure  and  saaple  porewater  pressure  responses  have  been 
plotted  together  on  each  figure.  The  'series'  designation  for  relative 
density  has  been  used  to  group  together  data  having  approxiaately  the 
saae  relative  density.  The  designations  include  data  that  is  within  10 
percent  greater  than  the  series  nuaber  (including  the  series  nuaber). 
For  exaaple,  a  '40%'  series  designation  would  include  all  data  for  a 
relative  density  fron  40  percent  to  49  percent. 

The  pressure-tiae  histories  are  indicative  of  the  systea  response 
during  and  after  loading.  On  each  figure,  the  traces  of  the  confining 
pressure  and  saaple  porewater  pressure  follow  each  other  closely  in 
their  response  trends.  The  two  curves  are  slightly  offset  froa  one 
another  in  the  tine  doaain  due  to  the  relative  locations  of  each 
transducer.  The  saaple  peak  porewater  pressure  is  greater  than  the 
applied  stress  peak  values  for  each  iapact. 

In  all  cases  the  confining  pressure  transducer  response  returned  to 
its  original  baseline  value  once  the  conpressive  stress  wave  energy  had 
dissipated.  The  response  of  the  confining  pressure  transducer  was  as 
expected  since  the  static  confining  pressure  should  be  constant  if  the 
systea  is  not  allowed  to  drain.  The  residual  excess  porewater  pressure 
indicated  by  the  saaple  transducer,  was  above  its  original  baseline 
value  after  each  loading  and  continued  to  increase  with  each  successive 
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impact.  The  response  of  the  sample  porewater  transducer  was  as 
expected  since  an  increase  in  the  residual  excess  porewater  pressure 
should  be  maintained  for  undrained  conditions.  Liquefaction  occurs 
when  the  residual  exoess  porewater  equals  the  effective  stress  which  is 
elso  when  the  back  pressure  plus  the  residual  excess  porewater  pressure 
equals  the  confining  pressure. 

The  information  from  the  pressure-time  history  records  was  used  in 
performing  the  data  analysis  presented  in  Chapter  VI.  The  porewater 
pressure  response  was  evaluated  as  a  function  of  initial  effective 
stress,  initial  sample  density,  peak  compressive  strain. 


j  D.  Composite  Data  Curves 

jj  Figures  5.29  through  5.33  show  the  porewater  pressure  ratio  as  a 

function  of  the  sum  of  the  sample  peak  porewater  pressures  grouped  by 
initial  relative  density.  Figure  5.34  shows  the  results  combined  on  a 
single  plot.  The  general  trend  is  for  the  porewater  pressure  ratio  to 
decrease  with  increasing  initial  relative  density  at  a  given  sample 
peak  porewater  pressure. 

Similarly,  Figures  5.35  through  5.38  show  the  porewater  pressure 
ratio  as  a  function  of  the  sum  of  the  peak  sample  porewater  pressures 
grouped  by  initial  effective  stress.  Figure  5.39  shows  the  results 
combined  together  on  a  single  plot.  Here,  the  general  trend  is  for  the 
porewater  pressure  ratio  to  decrease  with  incressing  initial  effective 
stress  at  a  given  sample  peak  porewater  pressure. 

It  is  evident  from  Figures  5.29  through  5.39  that  both  the  initial 
relative  density  and  initial  effective  stress  sre  important  factors 
that  influence  porewater  pressure  increases  in  a  saturated, 
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cohesionless  soil.  The  analysis  presented  in  Chapter  VI  develops 
empirieal  relationships  thst  relate  the  porevater  pressure  ratio  to 
these  two  parameters  and  the  peak  compressive  strain. 


Figure  5.9  Pressure-Time  Histories  for  D  =  *40%'  Series  and  o'  =  172  KPa  (Low  Impact 
Stress-First  Impact) 


Pressure-Time  Histories  for  D  =  ’40%’  Series  and  c'  =  172  KPa  (Low  Impact 
Stress-Second  Impact)  r  ° 
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VI.  ANALYSIS  AND  APPLICATION  OF  RESULTS 


A.  Bulk  Modulus  and  Compressive  Stre»»  Wave  Propagation 

The  water  saturated  sand  samples  used  in  this  experimental 
investigation  represent  a  two-phase  medium  for  stress  wave  propagation. 
The  presence  of  the  solid  sand  particles  in  the  water  increases  the 
fluid  density  and  also  affects  the  compressibility  of  the  soil-water 
mixture.  The  changes  in  fluid  density  and  mixture  compressibility  need 
to  be  considered  in  evaluating  the  compressive  stress  wave  propagation 
velocity  in  a  mixture.  A  procedure  presented  by  Richart  et  al.  (1970) 
was  followed  and  will  be  outlined  here. 

The  total  mass  density,  pt  ,  of  a  aoil-water  mixture  can  be 
determined  from: 


P 


t 


pw 


+  e 
+  e 


(Eq.  6.1) 


where,  G#  is  the  specific  gravity  of  the  solid  particles,  pw  is  the 
mass  density  of  water  and  e  is  the  void  ratio. 

The  compressibility  of  a  soil-water  mixture,  considering  the  solid 
particles  to  be  suspended  in  the  water,  consists  of  two  contributing 
factors:  the  compressibility  of  the  solid  particles  and  the 
compreasibil ity  of  the  fluid.  The  mixture  compressibility  can  be 


determined  from: 
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B 


mix 


B 


‘1  +  e' 


B 


1  +  •' 


(Eq.  6.2) 


where,  B^  it  the  bulk  nodules  of  the  weter  end  B#  is  the  bulk  atodulus 
of  the  soil  particles.  For  quartz  particles,  Bg  is  about  30680  MPa  and 
for  distilled,  de-aired  fresh  water  at  20  degrees  Celsius,  B^  is  about 
2140  MPa  (Richart  et  al.,1970).  It  should  be  noted  that  the  theory  of 
mixtures  assumes  that  the  solid  particles  are  suspended  in  the  water. 
However,  the  total  stress  will  be  soaiewhat  larger  since  the  solid 
particles  are  actually  in  contact  with  one  another.  For  small  strain 
conditions,  this  difference  is  small  and  can  be  neglected. 

Using  the  bulk  modulus  and  the  total  mass  density  for  the  mixture, 
the  compressive  stress  wave  propagation  velocity  through  the  mixture, 
V  .  ,  can  be  found  from: 


mix 


(Eq.  6.3). 


The  value  of  V  ^  includes  the  effects  of  density  and  compressibility. 
Calculated  values  for  the  Monterey  No.  0/30  sand  used  in  this 
experimental  investigation  are  given  in  Table  6.1.  Figure  6.1  shows 
the  acoustic  impedance  as  a  function  of  void  ratio  and  Figure  6.2  shows 
the  compressive  stress  wave  propagation  velocity  as  a  function  of  void 
ratio. 


B.  Peak  Particle  Velocity  and  Peak  Strain 

Two  values  of  particular  interest  thst  are  useful  in  interpreting 
the  results  of  this  investigation,  are  the  peak  particle  velocity  and 
the  peak  compressive  strain.  For  saturated,  sandy  soils  the  peak 
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Table  6.1  Streaa  Vave  Propagation  Parameters  for  Monterey 
No.  0/30  Sand 


D 

r 

<%) 

e 

pt 

(Kg/M3) 

B  .  V  . 
mix  mix 

(KPa) 

V  .  /V 

mix  w 

(M/sec) 

>t  Vc 

(Kg/O^-sec) 

0 

.803 

1915 

4421407 

1519 

1.013 

2908885 

10 

.779 

1926 

4485912 

1526 

1.017 

2940602 

20 

.755 

1940 

4554171 

1532 

1.021 

2972080 

30 

.731 

1953 

4626521 

1539 

1.026 

3005667 

40 

.707 

1967 

4703343 

1546 

1.031 

3040982 

50 

.683 

1980 

4785064 

1555 

1.037 

3078900 

60 

.659 

1995 

4872167 

1563 

1.042 

3118185 

70 

.635 

2009 

4965202 

1572 

1.048 

3158148 

80 

.611 

2024 

5064797 

1582 

1.055 

3201968 

90 

.587 

2040 

5171672 

1592 

1.061 

3247680 

100 

.563 

2056 

5286654 

1604 

1.069 

3297824 

Note:  V  *  the  compressive  stress  vave  velocity  in  fresh  water 
(Vw  *  1500  M/sec  at  20  degrees  Celsius) 


Figure  6.1  Acoustic  Impedance  as  a  Function  of  Void  Ratio  for  Monterey  No.  0/30  Sand 
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particle  velocity.  ,  can  be  determined  by  using  Equation  6.3  and 

rewriting  Equation  2.19  in  term*  of  the  peak  porewater  pressure 

response,  u  .  ,  as  follows: 

P* 


V  ‘  <»*W  V 


<E«.  6.4) 


The  peak  particle  velooity,  Vpk  .  can  then  be  obtained  by  rearranging 
Equation  6.4: 


where.  is  calculated  from  Equation  6.3. 

The  peak  coaipressive  strain,  e^  ,  can  then  be  determined  by 
combining  Equations  2.19  and  6.3  as  follows: 


la  .  -.i 

T.  ,t  iV2 


(Eq.  6.6). 


Tables  6.2  through  6.12  summarize  the  results  of  this  experimental 


investigation.  In  general,  the  experimental  results  show  that 

liquefaction  occurred  when  the  sum  of  the  peak  compressive  strains 

exceeded  about  lxlO-1  percent.  When  the  sum  of  the  peak  compressive 

_2 

strains  was  less  than  about  1x10  percent,  substantial  porewater 


pressure  increases  did  not  occur.  These  values  suggest  threshold 


compressive  strain  limits  for  Monterey  No.  0/30  sand  at  the  relative 
densities  and  effective  stresses  examined  in  this  experimental 
investigation. 

It  was  possible  to  obtain  porewater  pressure  ratios  greater  than 
one  in  several  of  the  experimental  results.  Other  researohers  have 


137 

j 

Table  6.2 

Peek  Porevater  Pressures  end  Peek  Compressive  Strains  for  f 

Monterey  No.  0/30  Send  for  Low  Impact 

Stress  Loading  and  i 

• 

Dr  -  '0**  Series 

j 

m 

j 

Test  I.D 

o'  D 

o  r 

Iapact 

u 

pk 

.  <2) 

pk 

PPB(1) 

(XPs)  («) 

(XPa) 

(%) 

SOOX* 

86  10.0 

1 

272 

.00606 

.329 

2 

118 

.00263 

.457 

3 

235 

.00524 

.743 

4 

108 

.00241 

.735 

5 

335 

.00747 

1.078  | 

6 

145 

.00323 

1.134  j 

soon 

172  4.6 

1 

461 

.01030 

.774  | 

2 

217 

.00485 

.902  1 

3 

362 

.00808 

1.114  J 

4 

443 

.00990 

1 .220  J 

5 

643 

.01435 

1.215  Q 

6 

588 

.01314 

1.219  | 

S00X2 

345  7.5 

1 

1222 

.02734 

.557  1 

2 

588 

.01316 

.614 

.740 

3 

670 

.01498 

» 

4 

244 

.00546 

.742 

5 

516 

.01154 

.848 

6 

917 

.02052 

.956 

S00X4 

690  7.5 

1 

81 

.00182 

.020 

2 

543 

.01214 

.211 

3 

208 

.00466 

.225 

4 

670 

.01498 

.341 

5 

561 

.01256 

.378 

6 

2786 

.06232 

.873 

Note: 

(1)  Measured 

(2)  Calculated  fro*  Equation  5 
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Table  (.3  Pitk  Pormttr  Pr«<ior*i  and  Peak  Coapressive  Straiaa  for 
Moaterey  No.  0/30  Sand  for  Low  Iapact  Stress  Loading  and 
Dr  -  *20%'  Series 


Teat  I.D. 

o' 

o 

(KPa) 

D 

r 

<%> 

Iapact 

«vU) 

Pk 

(KPa) 

■  „<2> 

Pk 

<*> 

PPR(1) 

S20X* 

86 

29.2 

1 

353 

.00764 

.116 

2 

516 

.01116 

.467 

3 

316 

.00685 

.565 

4 

362 

.00783 

.635 

5 

634 

.01371 

.576 

6 

379 

.00821 

.634 

S20X1 

172 

27.9 

1 

524 

.01136 

.105 

2 

416 

.00903 

.423 

3 

452 

.00981 

.740 

4 

469 

.01017 

.743 

5 

226 

.00490 

.849 

6 

452 

.00981 

.889 

S20X2 

345 

29.6 

1 

661 

.01429 

.229 

2 

489 

.01057 

.117 

3 

1048 

.02267 

.132 

4 

570 

.01233 

.213 

5 

1276 

.02760 

.536 

6 

914 

.01977 

.609 

S20X4 

690 

28.8 

1 

380 

.00823 

.185 

2 

441 

.00956 

.181 

3 

425 

.00921 

.278 

4 

543 

.01175 

.384 

5 

516 

.01117 

.422 

6 

272 

.00588 

.423 

Note:  (1)  Measured 

(2)  Calculated  froa  Equation  5.6 
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T»bl»  6.4  Peak  Pornitir  Pressures  sad  Peak  Conpressive  Strains  for 
Monterey  No.  0/30  Sand  for  Low  Inpact  Stress  Loading  and 
Dr  -  *40%'  Series 


Teat  I.D. 

o* 

0 

(IPs) 

D 

r 

(%) 

Iapact 

u 

pk 

(IPs) 

•  <2> 
*pk 

(%) 

ppRU) 

S40X* 

86 

46.7 

1 

517 

.01087 

.421 

2 

398 

.00836 

.567 

3 

448 

.00942 

.527 

4 

769 

.01617 

.898 

5 

706 

.01484 

1.064 

6 

996 

.02093 

1.222 

S40X1 

172 

47.1 

1 

1692 

.03554 

.606 

2 

353 

.00742 

.656 

3 

416 

.00875 

.695 

4 

1140 

.02396 

.881 

5 

407 

.00856 

.935 

6 

443 

.00931 

.996 

S40X2 

345 

45.9 

1 

54  3 

.01142 

.062 

2 

480 

.01010 

.114 

3 

3061 

.06444 

.531 

4 

751 

.01581 

.562 

5 

878 

.01848 

.572 

6 

2006 

.04224 

.691 

S40X4 

690 

46.7 

1 

724 

.01522 

.345 

2 

136 

.00286 

.347 

3 

597 

.01255 

.447 

4 

570 

.01199 

.488 

5 

760 

.01597 

.527 

6 

525 

.01103 

.551 

Note:  (1)  Measured 

(2)  Calculated  fron  Equation  5.6 
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Table  6.5  Peek  Poreuater  Pressures  and  Peak  Coapressive  Strains  for 

Monterey  No.  0/50  Sand  for  Lou  Inpact  Stress  Loading  and 

D  -  *60%'  Series  (First  Set) 
r 


Test  I.D. 

o’ 

0 

(XPa) 

D 

r 

<%) 

Inpact 

u  “> 
pk 

(XPa) 

«  <2) 

Pk 

(%> 

PPR(l> 

S60X* 

86 

No  Data 

Retrieved 

S60X1 

172 

66.7 

1 

941 

.01908 

.164 

2 

3592 

.07281 

.336 

3 

796 

.01614 

.403 

4 

796 

.01614 

.457 

5 

896 

.01815 

.525 

6 

480 

.00973 

.495 

S60X2 

345 

67.5 

1 

652 

.01319 

.216 

2 

823 

.01666 

.346 

3 

335 

.00678 

.346 

4 

832 

.01684 

.409 

5 

751 

.01520 

.399 

6 

4220 

.08540 

.560 

S60X4 

690 

65.4 

1 

543 

.01103 

.335 

2 

965 

.01961 

.369 

3 

534 

.01084 

.450 

4 

263 

.00534 

.439 

5 

697 

.01416 

.495 

6 

1339 

.02721 

.610 

Note:  (1)  Measured 

(2)  Calculated  from  Equation  5.6 
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Table  6.6 


Peak  Porevater  Pressures  and  Peak  Conpressive  Strains  for 
Monterey  No.  0/30  Sand  for  Low  Znpaot  Stress  Loading  and 
“  *60%'  Series  (Second  Set) 


Test  I.D. 

o' 

0 

(XPa) 

D 

r 

<%> 

Iapact 

»  “> 
pk 

(EPa) 

•  h<2> 

Pk 

(%) 

PPR(1) 

S60X+ 

86 

67.5 

1 

349 

.00767 

.464 

2 

190 

.00385 

.606 

3 

63 

.00128 

.562 

4 

308 

.00622 

.747 

5 

352 

.00712 

.741 

6 

398 

.00805 

1.021 

S60X1 

172 

67.1 

1 

380 

.00769 

.165 

2 

244 

.00494 

.263 

3 

498 

.01008 

.370 

4 

724 

.01466 

.411 

5 

652 

.01320 

.478 

6 

643 

.01302 

.527 

S60X2 

343 

66.3 

1 

290 

.00587 

.117 

2 

724 

.01469 

.298 

3 

416 

.00845 

.357 

4 

507 

.01028 

.345 

5 

950 

.01927 

.397 

6 

1077 

.02185 

.504 

S60X4 

690 

65.8 

1 

480 

.00974 

.185 

2 

325 

.00661 

.208 

3 

434 

.00882 

.225 

4 

1149 

.02333 

.354 

5 

1031 

.02094 

.410 

6 

No  Data  Recorded  this 

Impact 
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Table  6.7  Peek  Porewater  Pressure*  sad  Peek  Conpressive  Strains  for 
Monterey  No.  0/30  Send  for  Low  Impact  Stress  Loading  and 
Dr  -  *80%'  Series 


Test  I.D. 

o' 

0 

(IPa) 

D 

r 

(*) 

Impact 

u  <*> 
Pk 

(KPs) 

ak<2> 

Pk 

<%) 

PPR(1) 

S80X* 

86 

85.8 

1 

796 

.01554 

.636 

2 

914 

.01784 

.890 

3 

561 

.01095 

.932 

4 

1240 

.02418 

.960 

5 

1575 

.03072 

1.060 

6 

1057 

.02066 

1.073 

S80X1 

172 

87.9 

1 

950 

.01845 

.532 

2 

805 

.01564 

.619 

3 

2151 

.04178 

.798 

4 

579 

.01125 

.850 

5 

878 

.01705 

.902 

6 

878 

.01705 

.938 

S80X2 

345 

85.4 

X 

697 

.01361 

.364 

2 

715 

.01396 

.424 

3 

914 

.01785 

.468 

4 

643 

.01255 

.476 

5 

733 

.01431 

.532 

6 

661 

.01290 

.511 

S80X4 

690 

86.3 

1 

715 

.01394 

.331 

2 

416 

.00812 

.358 

3 

597 

.01164 

.399 

4 

389 

.00758 

.408 

5 

470 

.00917 

.435 

6 

461 

.00899 

.411 

Note:  (1)  Measured 

(2)  Calculated  from  Equation  S.6 
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Table  6.8  Peak  Porevater  Presaores  and  Peak  Conpreaaive  Strains  for 
Monterey  No.  0/30  Sand  for  High  Inpact  Stress  Loading  and 
D_  -  *0%'  Series 


Test  I.D. 

o' 

0 

(KPa) 

Dr 

(%) 

Inpact 

u  (1) 
Pk 

(KPa) 

.  <2> 
Pk 

(%) 

PPR(1) 

SOOX1 

172 

4.6 

1 

5711 

.12832 

.662 

2 

6674 

.14997 

.664 

3 

6730 

.15121 

.664 

4 

6219 

.13974 

.600 

S00X2 

343 

3.8 

1 

6585 

.14813 

1.077 

2 

6757 

.15201 

1.593 

3 

6985 

.15713 

1.276 

4 

5400 

.12148 

1.310 

S00I4 

690 

3.8 

1 

2912 

.06532 

1.007 

2 

2601 

.05835 

.994 

3 

4693 

.10527 

.994 

4 

4807 

.10780 

1.068 

Note:  (1)  Measured 

(2)  Calculated  fron  Equation  3.6 


-  No  loadings  were  done  for  a £  **  86  KPa 
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Table  6.9  Peak  Porevater  Preasorea  and  Peak  Compressive  Strains  for 
Monterey  No.  0/30  Sand  for  High  Impact  Stress  Loading  and 
Dr  =  »20%'  Series 


Test  I.D. 

o' 

0 

Dy  Impact 

(1) 

Upk 

.  <2> 

Pk 

PPR 

(IPs) 

(%) 

(EPa) 

(%) 

S20X1 

172 

23.8 

1 

3817 

.08332 

1.015 

2 

4580 

.09999 

.994 

3 

5033 

.10986 

.994 

4 

4750 

.10369 

.994 

S20X2 

345 

27.5 

1 

5400 

.11719 

.994 

2 

5655 

.12271 

.994 

3 

5768 

.12517 

.828 

4 

3421 

.07425 

.837 

S20X4 

690 

22.1 

1 

5768 

.12624 

1.035 

2 

6274 

.13733 

1.035 

3 

6247 

.13673 

.994 

4 

5089 

.11139 

.994 

Note:  (1)  Measured 

(2)  Calculated  from  Equation  5.6 


No  loadings  were  done  for  a £  »  86  KPa 


Table  6.10  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 

Monterey  No.  0/30  Sand  for  High  Impact  Stress  Loading  and 

D  -  '40%'  Series 
r 


Test  I.D. 

o' 

o 

(IPa) 

D 

r 

(%) 

Impact 

»  «> 
pk 

(KPa) 

a  <2> 

Pk 

(%) 

ppR(D 

S40X1 

172 

44.2 

1 

4524 

.09550 

.504 

2 

6647 

.14032 

.829 

3 

6902 

.14571 

.829 

4 

6619 

.13974 

.871 

S40X2 

345 

46.7 

1 

4693 

.09866 

.834 

2 

3676 

.07727 

.919 

3 

4015 

.08440 

.929 

4 

4269 

.08974 

.928 

S40X4 

690 

46.7 

1 

3845 

.08083 

.787 

2 

5570 

.11708 

.911 

3 

4976 

.10460 

.952 

4 

7212 

.15160 

1.076 

Note:  (1)  Measnred 

(2)  Calculated  from  Equation  5.6 


No  loadings  were  done  for  o^ 


86  KPa 
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Table  (.11  Peek  Porewater  Pressures  and  Peak  Compressive  Strains  for 
Monterey  No.  0/30  Sand  for  High  Impact  Stress  Loading  and 
Dr  -  *60%'  Series 


Test  I.D. 

o  * 

0 

(KPa) 

D 

r 

<%> 

Impact 

"v(1> 

pk 

(KPa) 

•  t<J) 
pk 

(%) 

PPRU) 

S60X1 

172 

61.1 

1 

4156 

.08504 

.834 

2 

5711 

.11686 

1.166 

3 

3421 

.07000 

1.166 

4 

4326 

.08851 

1.143 

S60X2 

345 

64.2 

1 

6022 

.12265 

.662 

2 

5117 

.10422 

.745 

3 

5598 

.11401 

.909 

4 

6019 

.12259 

.994 

S60X4 

690 

62.5 

1 

5994 

.12245 

.624 

2 

5681 

.11607 

.704 

3 

5771 

.11790 

.787 

4 

5599 

.11438 

.787 

Note:  (1)  Measured 

(2)  Calculated  from  Equation  5.6 

-  No  loadings  were  done  for  ■  86  KPa 


Table  6.12  Peak  Porewater  Preasures  and  Peak  Compressive  Strains  for 
Monterey  No.  0/30  Sand  for  High  Inpact  Stress  Loading  and 
D  *  '80% '  Series 


Note:  (1)  Measured 

(2)  Calculated  from  Equation  5.6 

-  No  loadings  were  done  for  **  86  KPs 
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also  reported  porevater  pressure  ratios  greater  than  one  in  the  field 


(Kok, 1977>1978b;  Stnder  and  Kok  1980,  Fragaszy  et  al.,1983)  and  in  the 


laboratory  (Kuygrok  and  Van  der  Eogel.1980,  Van  der  Kogel  et  al.,1981). 


A  ratio  greater  than  one  suggests  that  the  porevater  pressure  has 


ezeeeded  the  total  stress  and  the  effeotive  stress  has  becoae  negative. 


Theoretically,  the  porevater  pressure  ratio  should  not  exceed  a  value 


of  one.  Part  of  the  increase  can  be  explained  by  the  limitations  of 


the  electronic  measurement  and  recording  systems  part  can  be  attributed 


to  random  experimental  errors.  Beyond  these  reasons,  it  is  not  clear 


why  the  porevater  pressure  ratio  should  have  exceeded  unity. 


rrrrrmi  »tt smsi  rr?  mi 


The  ii 


in  porevater  pressure  in  a  saturated  soil  occurs  as  a 


result  of  the  tendency  for  the  soil  skeleton  to  decrease  in  volume 


under  an  applied  load.  As  the  volume  decrease  takes  place,  the  pore 


spaces  reduce  in  size  and  the  porevater  pressure  increases.  The 


potential  for  soil  skeleton  volume  changes  is  usually  quantified  in 


terms  of  the  relative  density,  ,  vhich  relates  a  soil's  density  to 


the  maximum  and  minimum  values.  For  soils  having  similar  particle 


types,  shapes  and  grain  size  distributions,  the  relative  density  is  a 


useful  index  parameter  that  can  be  estimated  in  the  field  by  a  standard 


penetration  test  (AS1M  D1586).  Soils  having  relative  densities  belov 


63  percent  have  been  observed  to  have  a  high  liquefaction  potential  for 


both  explosive-induced  ground  motions  and  seismic-induced  ground 


motions.  Hovever,  vhen  the  grain  characteristics  mentioned  differ 


considerably,  Lee  and  Fitton  (1968)  observed  differing  resistances  to 


liquefaction  at  the  same  relative  density  for  cyclic  loading  of 


saturated  sands.  Therefore,  in  sneh  eases  the  relative  density  data 
would  not  provide  a  useful  indioator  for  estimating  liquefaction 
potential.  Ishihara  and  Vatanabe  (1976)  explained  the  variations  in 
liquefaction  potential  in  tens  of  a  soil's  'volume  decrease  potential 
( VDP) . ' 

Figure  6.3  shows  two  elements  of  equal  volumes  for  two  different 
sands  at  the  same  relative  density.  The  soil  in  diagram  a,  has  a  much 
smaller  potential  for  volume  decrease  than  does  the  soil  in  diagram  b, 
even  though  they  are  at  the  same  relative  density.  The  potential  for 
volume  change  can  be  observed  by  noting  the  change  in  void  ratio 
possible  for  each  soil  to  achieve  100  percent  relative  density  packing. 
Based  on  their  relative  densities,  both  soils  would  appear  to  have  the 
same  potential  for  liquefaction.  However,  their  void  ratio  data 
indicates  the  contrary.  The  larger  the  volume  reduction  possible,  the 
larger  the  potential  porewater  pressure  increase  and  potential  for 
liquefaction.  Ishihara  and  Vatanabe  (1976)  present  the  volume  decrease 
potential  concept  in  equation  form  as: 


VDP 


e  -  e 


min 


(1  -  V  <‘..x 


°min^ 


(Eq.  6.7) 


where,  e  is  the  sample  void  ratio,  e  .  is  the  minimum  void  ratio,  e 

min  max 

is  the  maximum  void  ratio  and  is  the  relative  density  as  a  decimal. 
The  term  (e _ _  -  e  .  )  accounts  for  variations  in  sand  types  at  similar 

■tz  ■till 

relative  densities.  Whenever  information  on  void  ratio  limits  for  a 


granular  soil  is  available,  it  would  be  preferable  to  use  Equation  6.7 
to  estimate  volume  change  potential.  However,  in  the  field  often  the 
only  available  density  information  is  the  relstive  density.  Two 


■od«l«.  on*  using  relative  density  and  one  nsing  the  volue  decrease 
potential  concept,  have  been  proposed  and  are  presented  in  the  next 
section. 


d.  Miiiiiitiiii  toumitt  Aiilnii 

In  condncting  this  experinental  investigation,  three  iaiportsnt 
factors  that  inflnenoe  liquefaction  were  examined.  They  were,  the 
initial  effective  stress,  the  initial  sanple  density  and  the  applied 
compressive  stress  or  coapr essive  strain.  To  evaluate  the  influence  of 
each  factor,  a  statistical  analysis  of  the  data  was  performed  to 
develop  a  model  for  predicting  the  porevater  pressure  ratio,  PPR.  The 
porewater  pressure  ratio  ,  is  a  nondimensionalized  factor  defined  as: 

PPE  ,  °S  ~  °9  .  AH  .  ^Z.  (Eq.  6.8) 

0  o  o 

where, 

a'x  -  the  residual  effective  stress  after  the 
passage  of  the  compressive  stress  wave, 

o'  ■  the  initial  effective  stress, 
o 

An  ■  the  change  in  porewater  pressure  and  equal 
to  the  change  in  effective  stress,  and 

■  the  residual  excess  porewater  pressure. 

In  nondimensionalixed  form,  the  changes  in  porewater  pressure  due  to 
loading  of  soils  having  differing  initial  effective  stresses  and 
relative  densities  can  be  compared.  The  results  from  laboratory  and 
field  investigations  can  also  be  compered. 


On*  aethod  of  ivilutisi  liquefaction  potential  aa  a  nuaerical 
result,  ia  to  exaaine  the  porewater  pressure  ratio  increase  frou  an 
applied  loading.  The  larger  the  ratio  becoues  the  closer  the  soil  is 
to  being  liquefied.  When  the  soil  has  liquefied*  the  porewater  pressure 
ratio  equals  unity.  The  ratio  oan  be  related  to  the  initial  effective 
stress*  initial  saaple  density  and  the  saaple  strain  frou  the  applied 
eoapressive  loading. 

Considering  these  factors,  a  series  of  predictor  awdels  was 
investigated  to  find  one  that  would  be  intuitively  correct  fron  an 
engineering  standpoint  and  would  also  fit  the  data  well  statistically. 
TVo  aodels  were  developed,  one  including  relative  density  and  one 
including  the  volwae  decrease  potential  which  would  provide  flexibility 
in  applying  the  aodels  to  other  laboratory  and  field  studies.  A 
linear,  aultivariate,  regression  an nlyai$  was  performed  using  a  nodal 
of  the  fora: 

log  (PPR)  -  log  a  +  b  log  c  +  d  log  f  +  g  log  h  (Eq*  6,9) 

where. 

PPR  -  the  porewater  ratio,  defined  as  the  ratio  of  the 
residual  excess  porewater  pressure  to  the  initial 
effective  stress, 

a  "  the  constant  coefficient, 

c  “  the  initial  effective  stress  in  KPa, 

f  ■  the  initial  relative  density  in  percent, 
or  the  voluae  decrease  potential, 

h  “  the  sua  of  the  applied  peak  eoapressive  strains 
in  percent,  and 

b,  d,  g  ■  the  exponents  for  each  nodel. 
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Tht  Umax,  multivariate,  regreaaion  analysis  vas  performed  using 

the  program  'MINITAB'  (Ryan  et  al. ,1976,1982) .  Table  6.13  presents  the 

six  models  used  in  analysing  the  data  and  Table  6.14  gives  tvo  selected 

statistieal  parameters  used  to  evaluate  the  models.  The  tvo  parameters 

are  the  standard  error  of  estimate,  S,  and  the  coefficient  of 

determination,  R  .  The  standard  error  of  estimate  is  a  measure  of  how 

much  the  measured  value  deviates  from  the  average  value  predicted  by 

the  regression  equation.  In  other  words,  S  is  a  measure  of  the 

accuracy  of  the  predictions.  The  nailer  the  deviation  the  closer  the 

equation  predicts  the  data.  If  the  standard  error  of  estimate  was 

2 

xero,  then  the  equation  would  match  the  data  exactly.  The  R  value  is 

a  measure  of  how  much  of  the  variation  in  the  data  is  explained  by  the 

2 

regression  equation.  The  higher  the  R  value,  the  closer  the 
prediction  of  observed  behavior.  If  the  coefficient  of  determination 
vas  100  peroent,  then  all  of  the  data  could  be  explained  by  the 
regression  equation.  Based  on  the  results  of  the  statistical  analysis 


presented  in  Table  6.14,  tvo  models  were  selected.  In  terms  of  the 
volume  decrease  potential,  the  porevater  pressure  ratio  predictor  model 
has  the  following  form: 


PPR  -  (11.39)  (Zs  v)*321  (o')”*305  (VDP)*149  <Eq-  6'10) 

px  o 


where,  is  in  percent  and  o'  is  in  KPa.  In  terms  of  relative 
px  o 

density,  the  porevater  pressure  ratio  predictor  model  has  the  following 
form: 


Vr  - 


Table  <.13  Experimental  Porevater  Pres cure  Ratio  Predictor  Models 


Model 


(c)  PPM  -  (9.41)  (r  e.)*315  (H?)'*238 

P*  °Q 

(d)  PPM  -  (12.71)  (E  a  v)*325  (o'  D  )"*232 

pk  or 

(e)  PPM  -  (11.39)  (Z  a  .)*321  (o')"*305  (VDP)*149 

pt  o 

(£)  PPR  -  (16 .30)  (la.)’331  (o')"*308  (D  )"*179 

px  o  r 

Note:  -  a  .  and  D_  are  in  percent 
P*  1 

-  o'  is  in  KPa 
o 
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Table  6.14  Coeffielent  of  Determination  and  Standard  Error 
of  Batiaate  for  Porewater  Pressure  Ratio 
Predictor  Models  in  Table  6.13 


Model 

R2  (*) 

S 

a 

55.0 

.206 

b 

61.1 

.193 

0 

57.6 

.207 

d 

63.8 

.191 

e 

59.8 

.205 

f 

65.9 

.187 

Note:  R*  -  the  coefficient  of  determination 
S  -  the  standard  error  of  estimate 
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PPK 


(16.30)  (Zs  .)*331  (a*)"*308  (D  )'*179  lBq '  6A1) 

pi  o  r 


where,  D  end  •  .  ere  both  in  percent,  end  o'  it  in  EPi, 
r  pi  o 

The  equations  presented  in  Table  6.13  are  the  statistical  best  fits 
for  all  the  experimental  data  combined.  The  term  represents  the 
total  strain  applied  to  the  sample  at  a  given  impact  loading  including 
the  strain  from  previous  impacts.  Therefore,  the  Ii  .  quantity  is  a 
cumulative  strain  value  for  a  sample.  Equation  6.11  represents  the 
statistical  best  fit  of  all  the  models  examined  and  has  been  selected 
for  use  in  the  analysis  presented  in  this  chapter.  Figure  6.4  shows 
the  variation  of  residuals  about  the  regression  line  for  Equation  6.11. 
Tablea  6.15  and  6.16  give  compressive  strains  predicted  by  Equation 
6.11  for  a  porewater  pressure  ratio  of  0.25  and  1.0  respectively. 

It  would  be  useful  if  the  emulation  of  measured  peak  compressive 
strains  could  be  used  to  represent  the  strain  from  a  series  of 
successive  inspect  loadings  or  a  single  large  impact  loading.  To 
investigate  this,  the  two  models  given  by  Equations  6.10  and  6.11  were 
evaluated  considering  only  the  data  from  the  first  impact  loadings  from 
each  data  set.  In  terms  of  the  volume  decrease  potential,  the 
porewater  pressure  ratio  predictor  model  has  the  following  form: 

PPR  -  (10.90)  (s  .)*321  (o')'*305  (VDP)*149  (Eq*  6*12) 
px  o 


where,  a  .  is  in  percent  and  o’  is  in  KPa.  In  terms  of  relative 
px  o 

density,  the  porewater  pressure  ratio  predictor  model  has  the  following 
form: 


PPK 


(16.00)  («pk)‘331  (o;)'*308  <Dr)"*179 


(Eq.  6.13) 


(SUM  PK  STRRIN)  (ES)  -  (Dr) 


Figure  6.4a.  Variation  of  Residuals  About  the  Regression  Line  for 
Equation  6.11. 


SUM  PK  STRRIN  (50 


Figure  6.4b.  Porewater  Pressure  Ratio  for  Equation  6.11  as  a  Function 


of  Density  and  Effective  Stress. 
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Tibi*  6. IS  CoapressiTe  Strains  Predicted  by  Equation  6.11  for 
PPB  ■  0.2S  (Monterey  No.  0/30  Sand) 


Dr  <%> 


c^  (IPs) 

10 

20 

40 

60 

80 

86 

.00072 

.00105 

.00153 

.01908 

.00223 

172 

.00138 

.00201 

.00292 

.00364 

.00425 

34S 

.00264 

.00384 

.00558 

.00695 

.00812 

690 

.00503 

.00731 

.01064 

.01325 

.01548 

Note:  Predicted  strains  are  in  percent 


where,  0(  end  s^  ere  both  in  percent,  end  o'o  it  in  KPa.  The  stenderd 
error  of  estiaete  for  Equation  6,12  vet  0.193  end  the  stenderd  error  of 
eetiaete  for  Equation  6.13  wee  0.186. 

In  eonpering  Eqnetions  6.10  end  6.11  with  6.12  end  6.13,  it  cen  be 
seen  thet  they  vary  slightly  in  their  coefficients  end  thet  the 
stenderd  deviations  ere  nesrly  the  sene.  This  leeds  to  the  conclusion 
thet  Eqnetions  6.10  end  6.11  ere  velid  for  use  with  the  conpressive 
etrein  froa  e  single  lerge  loeding  or  froa  e  succession  of  sasller 
loedings  lending  to  the  seae  totel  peek  conpressive  strein  velue. 

Intuitively,  these  four  Eqnetions,  6.10,  6.11,  6.12  end  6.13 
correctly  represent  the  observed  behevior  in  thet  the  poreweter 
pressure  retio  veried  inversely  with  the  initiel  effective  stress  end 
reletive  density,  end  directly  with  the  voluae  decreese  potent iel.  As 
wes  shown  in  Figures  5.4,  5.5,  5.6  end  5.7  end  Teble  5.2,  the  soil 
skeleton  stiffness  elso  veried  directly  with  incressing  initial 
effective  stress  and  reletive  density  end  inversely  with  voluae 
decrease  potential. 

2.  tta jail  fomp»*itl9B  of  Madfll 

The  initiel  etteapt  of  the  aodel  was  to  predict  the  poreweter 
pressure  retio  as  e  function  of  peek  conpressive  strain  end  reletive 
density  or  voluae  decrease  potential  for  each  date  set  at  e  given 
initiel  effective  stress.  The  results  froa  the  initiel  aodel s 
indicated  e  definite  trend  with  variations  in  the  initiel  effective 


stress.  However,  these  results  gave  four  separate  equations  of  the 
seae  fora  with  different  constants  end  were  only  applicable  to  the  four 
particular  initiel  effective  stresses  eseaiued.  Tlerefore,  snore 


general  font  of  the  model  was  developed  to  include  the  initial 
effective  stress  as  a  variable. 

The  variation  in  initial  effective  stress  and  initial  saatple 
density  can  be  vieved  as  an  extension  of  the  relationships  into  s  third 
distension  as  is  shown  in  Figure  6.5.  In  the  figure,  it  can  be  seen 
that  a  three-dimensional  curved  planar  surface  extends  along  the  third 
axis.  The  surface  does  not  have  constant  end  points  along  the  extended 
edges.  For  a  given  initial  effective  stress  and  initial  density,  the 
three-dimensional  surface  becomes  a  plane  whose  abscissa  is  the  peak 
compressive  strain  and  ordinate  is  the  porewater  pressure  ratio. 

E.  Finite  Difference  Approximation 

A  finite  difference  solution  approximation  was  investigated  to 
simulate  the  laboratory  measured  pressure-time  histories  using  the 
equations  for  a  one-dimensional,  confined  compressive  loading. 
Equations  were  developed  to  model  both  the  transient  and  long-term 
porewater  pressure  response  in  the  soil  mass.  These  equations  are 
based  on  the  developments  presented  by  Biot  (1956a, b, 1962)  and  a 
general  approach  given  by  Ishihara  (1967).  The  finite  difference 
recursion  formulas  used  were  adapted  from  Biggs  (1964)  for  a  single 
degree  of  freedom,  viscoelastic  system.  Appropriate  modifications  were 
made  to  the  solution  to  model  the  soil  as  a  two-phase  materiel  and  to 
include  the  experimental  boundary  conditions. 

i.  Madil 


The  rheologic  model  consists  of  masses,  springs  and  dashpots  used 
to  represent  the  applied  compressive  stress  wsve  propagation  through 


the  confining  pressure  tube,  the  saturated  soil  sample  and  the  PVC 
momentum  trap.  Figure  6.6  shows  a  schematic  of  the  mass-spring-dashpot 
system.  The  confining  tube  was  treated  as  a  constrained  'bar'  of  water 
having  three  mass  elements  connected  by  linear  elastic  springs  whose 
stiffness  corresponded  to  the  bulk  modulus  of  water.  The  PVC  momentum 
trap  was  treated  in  the  same  way  as  the  confining  tube  using  two 
elements  and  the  bulk  modulus  for  PVC. 

In  modeling  the  soil-water  system,  consideration  was  given  to 
including  the  stiffness  of  the  mixture,  the  stiffness  of  the  soil 
skeleton  and  material  damping.  The  bulk  modulus  of  the  mixture  was 
evaluated  as  described  in  section  A  of  this  chapter.  The  soil  skeleton 
was  treated  in  two  parts,  with  a  separate  value  for  loading  and 
unloading  to  account  for  permanent  deformations  of  the  soil  skeleton. 
The  skeleton  stiffness  values  were  determined  from  one-dimensional, 
confined  compression  loadings  of  dry  Monterey  No. 0/30  sand  as  described 
in  section  A,  part  5  of  Chapter  V.  A  viscous  damping  term  was  also 
included  to  account  for  the  compressive  stress  wave  energy  dissipation 
due  to  system  and  material  damping. 

2.  Finite  Difference  Solution  Formulation 

The  finite  difference  approximation  uses  Biot's  (1962)  theory  for 
stress  wave  propagation  in  a  porous  elastic  medium.  Biot's  theory  was 
coupled  with  a  soil  skeleton  model  that  decreases  in  volume  whenever 
the  applied  compressive  loading  produces  a  compressive  strain  that 
exceeds  a  threshold  value. 

Biot's  modified  general  equations  for  a  porous  medium  have  the 


following  forms  for  the  component  parts: 


F(t)  Water  Soil -Water  PVC  /  Energy 

Elements  Elements  Elements  /  Absorbing 

^  Boundary 
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a.  overall  system  equilibrium: 


°ij.J  +  pt  *i 


pt  tti  +  >w  wi 


(Eq. 


b.  liquid  equilibrium: 


-p.i  *  *1 


*1 


••  B  ■*  ®i  ^ i 

U  +  *  0  w  +  "  *■  — * 1 


u  *Sr  wi  '  k 


(Eq. 


c.  porous  solid  equilibrium: 


ij.J 


+ 


(Eq. 


d.  stress-strain  (fluid  and  porous  solid): 


#ij 


+ 


(Eq. 


e 


”  wi,i 


(Eq. 


where. 


«>w 

wi 


the  total  stress  tensor, 
the  total  mass  density, 

the  components  of  the  gravitational  constant, 
the  acceleration  vector  of  the  solid  skeleton, 
the  mass  density  of  the  fluid, 
the  velocity  vector  of  the  fluid. 


(.14) 


(.15) 


(.1() 


(.17) 


(.18) 


wi 

8 


k 


6i.J 


*ij 

u 

9 


tko  acceleration  rector  of  tko  fluid, 

•  kiaotio  energy  correct ion  factor, 

tko  porosity  of  tko  soil  skoloton, 

Darcy's  eooffieiont  of  permeability, 

tko  effective  atroaa  tensor, 

tko  porevater  pressnro  tensor, 
tko  Eroaeeker  delta, 

tko  strain  tensor  of  tko  solid  skoloton, 

tko  displacement  rector  of  tko  solid  skeleton, 

tko  variation  of  tko  fluid  oontont 

tko  displacement  vector  of  tko  water  vitk  respect  to 
tko  solid  skoloton 


Por  tko  ease  of  one-dimensional  strain  vitk  undrained  loading 
conditions.  Equations  6.14,  6.15  and  6.16  reduce  to; 

dfl  _  „  ~  (Eq.  6.19) 

ds  pt  * 


_  m  n  ~  (Eq.  6.20) 

dz  pw  * 


o'  -  o  +  P 


(Eq.  6.21). 


Ike  stress-strain  conditions  (Equations  6.17  and  6.18)  reduoo  to: 
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&o '  -  Z(D„  AeP)  (Eq.  6.23) 

c 

vkirt, 

Dq  «  the  constrained  aodnlns  for  loading  or 
unloading  of  the  toil  skeleton. 

B  .  -  the  bulk  aodnlns  of  the  soil-water  aixtnre 

calculated  froa  Equation  (.2. 

T 

As  “  the  wolnaetrie  strain  of  the  aix,  and 

is5  -  the  inelastic  wolnaetrie 
strain  of  the  aix. 

The  following  reeurrenee  foranlas  presented  by  Biggs  (1964).  were 
used  to  ewalnate  the  teras  for  diaplaeeaent,  weloeity  and  acceleration: 

a.  displaceaent  (y): 

y(s+l)  .  2y(s)  _  y(s-l)  +  y(s)  (At)2  (Eq.  6.24) 


(Eq.  6.25) 


e.  acceleration  (y): 


[fI0..-  li%l  -  C£ILr-(,),  -_y(*~1))j 
-C.)  .  1 _ “ _ 1 


M(,)  +  C  (At)2 


(Eq.  6.26) 


where, 

s  =  the  time  at  any  atap  in  the  iteration. 
At  *  the  tine  increment  for  the  iteration, 
K  *  the  equivalent  spring  stiffness, 

C  ■  the  visoons  damping  coefficient, 

F  “  the  applied  force,  and 

M  -  the  mass. 


The  viscous  damping  coefficient,  C,  was  calculated  for  the  system 
at  each  time  step  from: 


C  -  (2)  (D)  [(I)  (M)]0*5 


(Eq.  6.27). 


The  damping  ratio,  D,  was  evaluated  from  the  pressure-time  history  data 
of  the  laboratory  experiments  and  hed  an  average  value  of  0.0167. 

F.  Comparison  of  Solution  With  Experimental  Results 

The  finite  difference  approximation  was  used  to  compare  a 
theoretical  analysis  with  observed  behavior  in  the  laboratory.  A 
series  of  computer  runs  using  a  triangular  loading  pulse  was  performed 
for  the  first  impact  loadings  of  two  laboratory  samples  at  two 
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different  initial  effective  stresses.  The  initial  relative  densities 
need  were  40  and  80  percent  and  the  initial  effective  stresses  used 
were  172  KPa  and  (90  KPa.  The  respective  loading  and  unloading 
constrained  siodnli  were  obtained  from  Figures  5.4,  5.5,  5.(  and  5.7. 
The  nodnli  were  detemined  at  straina  corresponding  to  the  experimental 
laboratory  results  and  their  nuaerieal  values  are  given  in  Table  5.2. 

Figures  (.7,  (.8,  6.9,  and  (.10  show  the  results  of  the  finite 
difference  analysis.  In  comparing  these  figures  with  the  experimental 
pressure-time  histories  presented  in  section  D  of  Chapter  V,  similar 
behavior  is  evident  for  the  porewater  pressure  response  both  during  and 
after  the  applied  loading.  There  is  also  a  marked  increase  in  the 
residual  excess  porewater  pressure.  Table  (.17  presents  a  comparison 
of  the  finite  difference  approximation  results  and  the  experimental 
results.  Variations  in  the  results  oan  be  attributed  to  the  use  of  a 
static  constrained  modulus  for  the  soil  skeleton  ibn  the  model  and 
laboratory  boundary  conditions  that  could  not  be  exactly  simulated  by 
the  finite  difference  approximation. 

G.  Application  of  Results 

In  Chapter  II,  several  empirical  soaling  equations  were  presented 
for  predicting  the  peak  compressive  stress  developed  by  a  buried  charge 
as  a  function  of  radial  distance  from  the  detonation  point  and  charge 
weight.  Equations  (.4,  (.5  and  (.(  developed  the  relationship  between 
peak  compressive  strain  and  peak  compressive  stress.  Since  the  peak 
compressive  strain  can  be  estimated,  the  statistical  models  developed 
in  this  experimental  investigation  can  be  used  to  approximate  the 
resulting  porewater  pressure  ratio  as  a  function  of  effective  stress 
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Table  6.17  Comparison  of  Porewater  Pressure  Ratios  Predicted  by 

Finite  Difference  Approximation  and  Equation  (.11  with 
Experimental  Results  for  Monterey  No.  0/30  Sand 


Note:  (1)  Experimental  Reaulte 

(2)  Predicted  by  Finite  Difference  Analysis 

(3)  Predicted  by  Equation  (.11 
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ud  toil  density.  This  provides  s  way  to  estimate  the  extent  of 
regions  of  porewater  pressure  inereases  at  distanees  outward  frost  a 
bnried  charge.  An  exaaple  nsing  this  approach  follows. 


To  estiaate  the  peak  coapressive  strains  froa  a  bnried  oharge.  it 
was  necessary  to  nse  one  of  the  peak  particle  velocity  scaling  laws 
given  in  Chapter  II  with  Equation  6.6.  For  this  exaaple.  the  scaling 
law  selected  was  Equation  2.18  (Drake  and  Little.  1983).  When  Equation 
2.18  is  substituted  into  Equation  6.6.  the  following  equation  for  peak 
compressive  strain  is  obtained: 


*pk 


(3.6)  (f)  ( 


-2.35 


V 

c 


(Eq.  6.28) 


where  all  terns  are  as  previously  defined. 

If  a  charge  placeaent  depth  of  5  neters  is  used.  Equation  2.3  gives 

a  charge  weight  of  7  kilograas  to  insure  a  contained  explosion.  For 

this  depth  and  weight,  the  scaled  depth  of  burst  is  2.6  M/Kg and 

1/3 

since  this  is  greater  than  0.56  M/Kg  (Figure  2.3).  the  groundshock 
coupling  factor,  f.  be cones  1.0.  Assuming  a  saturated  soil,  the 
eoapressive  stress  wave  velocity,  VQ  ,  can  be  tsken  as  1500  neters  per 
second.  Substituting  these  values  into  Equation  6.28  gives: 


*pk 


(0.086)  (R) 


-2.35 


(Eq.  6.29) 


for  a  oonfined  7  kilograa  charge  in  a  saturated  soil 
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The  porewater  praaanra  ratio  predictor  model  chosen  for  this 
example  is  that  given  by  Equation  6.11.  Two  relative  densities  were 
nsed.  40  and  80  percent,  and  the  effective  atress  and  peak  compressive 
strains  were  varied  with  the  radial  distance.  R.  from  the  buried  charge 
depth.  When  the  strains  calculated  from  Equation  6.28  are  used  in 
Equation  6.11.  the  result  gives  the  porewater  pressure  ratio  as  a 
function  of  depth  and  horizontal  distance  from  the  charge. 

Since  the  porewater  pressure  ratio  variation  has  been  determined  at 
points  away  from  the  buried  charge  depth*  a  series  of  contours  can  be 
drawn.  These  contours  are  symmetric  about  the  charge  location. 
Figures  6.11  and  6.12  show  the  porewater  pressure  ratio  oontours  for 
the  two  relative  densities  examined.  These  results  estimate  the  tones 
of  potential  porewater  pressure  increase.  Figure  6.13  compares  the 
maximum  radius  of  liquefaction  near  the  ground  surface  as  predicted  by 
the  approach  presented  and  Equation  2.7.  The  differences  in  the 
estimated  maximum  radius  of  liquefaction  can  be  attributed  to 
variations  in  soil  type  and  boundary  conditions  between  the 
experimental  inveatigation  and  the  field  tests. 
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VII.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  jJmyUH 

The  cause  of  earthquake-induced  liquefaction  is  veil  understood  and 
docuaMnted.  Shear  stress  wave  propagation  is  the  primary  cause  of 
increases  in  porewster  pressures.  Although  considerable  research  has 
been  conducted  on  explosive-induced  ground  notions  and  effects,  little 
work  has  been  directed  towards  examining  the  phenomenon  of  explosive- 
induced  liquefaction,  hence,  this  phenomenon  is  not  well  understood. 
This  experiMntal  investigation  represents  an  effort  to  examine  and 
establish  a  fundamental  understanding  of  compressionally-induced 
liquefaction. 

A  new  experimental  soil  dynamics  facility  has  been  developed  to 
investigate  the  transient  and  long-term  porewater  pressure  response  of 
saturated  soils  as  a  function  of  initial  density,  initial  effective 
stress  and  applied  load  intensity.  The  facility  is  capable  of 
generating  oompressive  shock  loadings  on  the  order  of  35000  KPa  with 
millisecond  rise  times  to  peak  stress.  Ssmples  of  wster  saturated 
Monterey  No.  0/30  sand  were  prepared  at  relative  densities  of 
approxiawtely  0  .  20,  40,  60  and  80  percent  and  effective  stresses  of 

86,  172,  345  and  690  KPa.  Each  sample  was  subjected  to  one- 
dimensional,  confined,  coaipressive  loadings  under  undrained  conditions. 
Experimental  results  have  shown  that  it  is  possible  to  liquefy  Monterey 
No.  0/30  sand  under  these  conditions.  Significant  porewater  pressure 
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increase*  were  observed  even  st  high  densities  and  high  effective 
stresses. 

It  is  evident  frost  this  investigation  that  liquefaction  becomes 

■ore  difficult  to  induce  with  iucreasing  density  and  effective  stress, 

sinoe  greater  energy  is  required  for  porewater  pressure  increases.  The 

intensity  of  the  applied  load  is  iaportant  and  liquefaction  oould  be 

induced  by  a  single  high  aaplitude  compressive  loading  or  a  series  of 

less  intense  compressive  loadings.  Threshold  compressive  strain 

values,  below  which  substantial  porewater  pressures  should  not  occur, 

—2 

were  determined  to  be  less  than  about  1x10  percent.  Liquefaction  was 
generally  observed  at  compressive  strains  greater  than  about  lxlO**1 
percent. 

The  data  analysis  has  provided  empirical  models  that  can  be  used  to 
estisute  the  compressive  strain  required  to  produce  liquefaction  in  a 
saturated  sand  as  a  function  of  density  and  effective  stress.  One 
model.  Equation  6,11,  has  been  used  with  an  empirioal  field  explosive 
scaling  equation  to  develop  a  method  of  approximating  the  extent  of  the 
liquefied  zone  in  the  field  when  buried,  contained  charges  are 
detonated  in  a  saturated  soil.  The  results  can  be  applied  to  estimate 
the  locations  of  regions  having  various  porewater  pressure  increases  at 
points  away  from  the  explosion. 

A  finite  difference  approximation  was  used  to  model  the  porewater 
pressure  response  observed  in  the  experimental  investigation.  The 
approximation  considered  the  soil-water  mixture  as  a  two-phase  material 
and  included  the  nonlinear,  inelastio  behavior  of  the  soil  skeleton. 
The  results  were  generally  consistent  with  the  laboratory  behavior  and 
showed  a  marked  increase  in  the  residual  exoess  porewater  pressure. 
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Variations  in  the  raaalta  can  be  attributed  to  using  a  static 
constrained  node las  for  the  soil  skeleton  in  the  model  and  laboratory 
boundary  conditions  that  oonld  not  be  exactly  simulated  by  the  finite 
difference  approximation. 

A  number  ».  significant  faotors  that  influence  the  occurrence  of 
liquefaction  have  been  reviewed.  Both  field  and  laboratory 
investigations  using  explosives  have  provided  consistent  observations 
about  the  behavior  of  saturated  soils  subjected  to  transient  loadings. 
Empirical  scaling  equations  to  predict  peak  particle  velocities,  peak 
compressive  strains  and  peak  compressive  stresses  have  been  presented. 

B.  figtttlBiiqai 

Based  on  the  results  of  this  experimental  investigation  as  well  as 
the  discussion  and  review  of  other  work,  the  following  conclusions  were 
drawn.  These  conclusions  pertain  to  the  laboratory  investigation  of 
saturated  Monterey  No.  0/30  sand  under  the  experimental  conditions  as 
previously  defined  in  the  text  and  include  the  following. 

Liquefaction  can  be  induced  by  dynamic  one-dimensional,  confined 
compressive  loadings  under  undrained  conditions  at  the  relative 
densities  and  effective  stresses  used  in  this  investigation. 


The  experimental  results  show  that  it  is  possible  to  obtain 
significant  porewater  pressure  increases  even  at  high  densities 
and  effective  stresses. 


Liquefaction  can  be  induced  by  a  single  large  compressive  loading 
or  a  series  smaller  compressive  loadings. 


The  initial  density  and  initial  effective  stress  are  significant 
factors  that  influence  the  occurrence  of  liquefaction.  In 
general,  a  sample  at  a  high  initial  relative  density  and  high 
initial  effective  stress  requires  more  applied  energy  to 
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produce  liquefaction  than  doe a  one  at  a  lower  initial  effective 
street  and  initial  density. 


Substantial  increases  in  porevater  pressure  did  not  occur  at 
eoapressive  strains  less  than  about  lxlO**2  percent. 
Liquefaction  could  be  induoed  at  compressive  strains  greater 
than  about  1x10  percent.  Both  single  and  aultiple 

eoapressive  loadings  were  observed  to  produce  liquefaction. 


The  finite  difference  analysis  indicates  that  liquefaction  is 
dependent  on  the  unloading  constrained  nodulus  of  the  soil 
skeleton  which  is  in  turn  dependent  on  the  initial  density  and 
initial  effective  stress.  The  unloading  effects  for  the  soil- 
water  aixture  are  significant  and  liquefaction  appears  to  occur 
during  unloading  due  to  the  nonlinear,  inelastic  behavior  of 
the  soil  skeleton. 


The  statistical  awdels  presented,  when  used  with  an  empirical 
scaling  equation  for  field  explosive  events,  can  be  used  to 
approximate  the  zone  of  liquefaction  around  a  buried  charge. 
The  models  need  to  be  verified  and  modified  to  account  for 
specific  conditions  in  the  field. 


The  results  of  this  investigation  have  increased  the  state-of- 
the-art  in  understanding  compressionally  induced  liquefaction 
and  have  extended  the  data  base. 


In  considering  the  nature  of  explosions  in  the  field  and  the 
results  of  this  investigation,  the  following  conclusions  can  be  drawn 
relative  to  the  use  of  buried  charges  in  the  field. 

Current  engineering  designs  are  generally  concerned  with  the  peak 
particle  velocity  and  damage  produced  by  the  explosion  and  do 
not  consider  any  potential  porewater  pressure  increases  or 
liquefaction  effects. 


The  damage  from  an  explosion  in  the  field  may  be  disproportionate 
to  the  amount  of  energy  released  if  liquefaction  occurs.  The 
damage  would  then  include  the  direct  effects  of  the  explosion 
itself  and  liquefaction. 


Liquefaction  may  occur  at  distances  and/or  peak  particle 
velocities  less  than  those  associated  with  structural  damage. 
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An  estimate  of  field  liquefaction  xonation  can  be  Bade  by  using 
the  approach  outlined.  However,  snail  scale  field  explosive 
tests  are  reconaended  to  verify  and  appropriately  adjust  the 
equations  to  aore  closely  aatch  site  specific  conditions  in  the 
field. 


C.  Rocommendat ions 

The  results  of  this  investigation  and  the  infomation  currently 
available  in  the  literature  have  provided  insight  into  research  needs 
to  better  understand  the  phenoaenon  of  liquefaction  as  produced  by 
intense  coapressive  loadings.  Based  on  this,  a  number  of 
recosimendations  can  be  made. 


Research  efforts  should  be  extended  to  investigate  the  effects  of 
different  soil  types  and  varying  degrees  of  saturation  on 
liquefaction  potential  of  soils  under  compressional  loading. 
Attention  also  should  be  given  to  evaluating  both  two  and  three 
dimensional  loading  phenomenon. 


Develop  a  series  of  fully  instrumented,  systematically  controlled, 
saall  scale,  field  explosive  tests  under  controlled  conditions. 
This  should  involve  the  variation  of  initial  density,  effective 
stress,  saturation  and  type  of  loading.  Correlations  can  then 
be  aade  between  field  and  laboratory  results. 


The  statistical  stodels  presented  and  the  empirical  equations  for 
explosively  produced  peak  particle  velocities  should  be 
evaluated  with  respect  to  site  specific  conditions  and  modified 
accordingly  to  account  for  field  conditions. 


The  effects  of  liquefaction  and  soil  property  changes  should  be 
further  investigated  and  incorporated  into  engineering  analysis 
and  design  procedures  that  consider  explosive  loading  of  soils 
and  structures. 


Develop  finite  element  models  that  will  account  for  the  loading 
conditions  and  boundary  conditions  of  this  investigation  and 
extend  them  to  two  and  three  dimensional  analyses.  Since  many 
models  treat  the  soil-water  system  as  a  one-phase  material,  it 
will  be  necessary  to  evaluate  a  given  element  of  soil  as  a 


two-phase  and  three-phase  naterial  to  store  accurately  evaluate 
the  behavior  under  loading. 

Develop  siethods  for  evaluating  the  dynastic  constrained  aodulus  of 
a  soil  for  both  loading  and  unloading  to  be  used  in  nusterical 
aodels. 


Establish  a  standard  data  set  to  be  obtained  frost  all  field  tests 
designed  to  investigate  explosively-induced  liquefaction.  This 
should  include  a  complete  and  appropriate  set  of  information 
including  soil  index  properties,  sub-surface  profile 

characterization,  groundwater  table  location,  degree  of 
saturation  of  the  soil  and  local  geologic  conditions. 

Detailed  information  should  also  be  obtained  for  the  explosive 
charge  being  used  including  type,  location,  weight,  delay  times 
and  charge  patterns. 


Examine  available  field  ground  motion  records  from  explosive 
events  to  define  where  the  transition  from  one-dimensional 
loading  to  shear  loading  occurs  and  evaluate  these  findings  in 
terms  of  their  interaction  in  producing  liquefaction. 
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